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Chapter  I  .  j- 

INTRODUCTION  | 

H 

I  . 

Interlaminar  stresses  play  an  important  role  in  the  load  trans-  i, 

fer  mechanism  in  composite  laminates.  Both  numerical  and  experimental 
results  have  demonstrated  that  when  a  thin  laminate  is  subjected  to  a  j 

uniaxial  extension  (Fig.  1),  there  exist  highly  localized  stress  con¬ 
centration  regions  near  the  free  edges,  the  so-called  boundary  layer 
regions.  This  phenomenon  has  been  suggested  to  be  the  dominant  factor 
initiating  failure  of  some  composite  laminates. 

The  present  thesis  will  analyze  the  boundary  layer  by  "perturb¬ 
ing"  the  exact  elasticity  equations  with  a  stretching  transformation. 

Solutions  to  these  transformed  equations  provide  a  higher  order  analysis 

than  idealized  lamination  theory  [1].*  Hence  better  insight  into  the 

interlaminar  stress  behavior  is  obtained  using  the  perturbation 

analysis  [2].  i 

1.1  REVIEW  OF  LITERATURE  j 

Bogy  [3]  analyzed  a  bonded  material  containing  two  mutually  dis-  I 

h 

similar  orthogonal  wedges  under  arbitrary  tractions.  The  stress  fields  j 

were  found  to  contain  a  mathematical  singularity  at  the  intersection  of  | 

the  interfacial  plane  and  the  loaded  surface.  Hein  [4]  studied  the  j 

J 

residual  stresses  in  a  two-material  wedge  and  found  similar  behavior.  i 

I 

♦Numbers  in  brackets  jefer  to  the  references  l.isted  in  the  j 

bibliography.  j 
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Hess  [5]  developed  a  plane  elasticity  solution  for  the  end  problem  in 
a  two-layer  laminated  strip  and  showed  a  mathematical  singularity, 
defined  to  be  a  point  where  the  convergence  of  an  eigenfunction  expan¬ 
sion  could  not  be  attained.  Puppo  and  Evensen  [6]  modeled  the  finite- 
width  symmetric  laminate  as  a  set  of  anisotropic  layers  separated  by 
isotropic  shear  layers.  Each  anisotropic  layer  was  assumed  to  be  under 
generalized  plane  stress,  i.e.,  the  out-of-plane  normal  stress 
(Fig.  1)  is  zero  and  the  in-plane  stresses  and  displacements  are  the 
thickness  averages  of  the  actual  values.  Solutions  to  the  corresponding 
equilibrium  equations  showed  tiiat  while  the  interlaminar  shear  stresses 
vanish  everywhere  for  a  laminate  of  infinite  width,  they  attained 
maximum  finite  values  near  the  free  edge  of  a  finite  width  laminate. 
Furthermore,  in  regions  far  av/ay  from  the  free  edge,  the  solution 
agreed  well  with  the  classical  lamination  theory  flj*  A  complete  three 
dimensional  analysis  was  carried  out  by  Pipes  and  Pagano  [7]  using  the 
finite  difference  technique  to  solve  the  exact  elasticity  equations. 

The  results  showed  good  agreements  with  those  of  Reference  [6]  except 
at  the  free  edge  where  the  interlaminar  shear  stress  seemed  to 
grow  without  bound  for  some  laminates.  Due  to  the  approximate  nature 
of  the  finite  difference  analysis,  however,  no  evidence  was  available 
to  show  the  intensity  of  the  suggested  singularity.  Isakson  and  Levy 
[8]  used  a  finite  element  approach  to  analyze  a  model  similar  to  that 
of  Reference  f6].  Based  on  the  constant  strain  assumption  within  each 
element,  the  corresponding  stresses  were  obtained  from  the  constitu¬ 
tive  equations.  The  total  elastic  strain  energy  was  calculated  and 
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minimized  [9]  to  yield  a  set  of  simultaneous  linear  algebraic  equa¬ 
tions.  Levy,  et  al.  [10]  used  the  same  model  and  formulation  as 
Reference  [8]  to  further  investigate  the  elastic  and  plastic  inter¬ 
laminar  shear  deformations  in  the  laminate.  The  out-of-plane  "peel 
stress"  was  not  taken  into  account  in  both  studies  due  to  the  modeling. 
Results  from  these  solutions  were  quite  similar  to  those  of  Reference 
[6]  except  at  the  free  edge  where  the  interlaminar  shear  stress  tx2 
was  40%  lower  than  that  of  Reference  [6].  This  presumably  v/as  due  to 
the  limitations  of  the  finite  element  approximation.  Improvements 
were  made  by  Rybicki  [11]  who  carried  out  a  three-dimensional  finite 
element  analysis  based  on  a  complimentary  energy  formulation  in  terms 
of  three  Maxwell  stress  functions.  These  functions  resulted  in  a  set 
of  simultaneous  linear  algebraic  equations  which  were  solved  by  Gauss 
reduction  and  the  back  substitution  process.  The  "peel  stress"  was 
obtained  in  this  investigation.  The  results  showed  excellent  agreement 
with  References  [6]  and  [7]  in  regions  removed  from  the  free  edges, 
while  near  the  edges  the  interlaminar  shear  stress  ag>'eed  only  with 
Reference  [6];  the  magnitude  of  was  much  lower  than  the  singular 
value  of  Reference  [7].  The  approximate  nature  of  the  finite  element 
foraulation  for  the  laminated  plate  apparently  leads  to  questionable 
and  quite  possibly  poor  results  at  the  exact  free  edge.  Pipes  [12] 
used  the  finite  difference  procedure  to  carry  out  extensive  parametric 
studies  including  laminate  geometry,  fiber  orientations  and  stacking 
sequerces.  The  program  used  in  Reference  [12]  is  capable  of  handling 
no  more  than  an  8  layer  symmetric  laminate  owing  to  the  limited 
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computer  capacity. 

Several  attempts  have  recently  been  made  to  verify  experimentally 
the  numerical  predictions.  Results  by  Pipes  and  Daniel  [13], 

Herakovich  [14],  and  Oplinger,  et  al.  [15]  have  all  shov/ed  significant 
stress  ctncentration  behavior  near  the  free  edges.  Although  stress 
intensities  were  not  determined  in  these  studies,  there  were  strong 
evidences  to  support  the  numerical  prediction  of  significant  stress 
concentrations  near  the  free  edge. 

Pipes  and  Pagano  [16]  more  recently  developed  an  analytical 
solution  to  the  elasticity  equations  under  the  assumptions  of  zero 
interlaminar  noraal  stress,  o^,  and  zero  transverse  normal  stress,  oy, 
for  the  [±45]j  laminate.  Pagano  [17]  obtained  yet  another  approximate 
solution  following  the  cylindrical  bending  theory  of  Whitney  and 
Sun  [18].  Good  agreement  with  the  elasticity  solution  of  Reference 
[7]  was  found  for  the  interlaminar  normal  stress,  a^,  (the  "peel 
stress")  on  the  midplane  of  a  bidirectional  [O/QO]^  laminate.  However, 
the  solution  did  not  recognize  the  stress  free  boundary  conditions 
Tyz(±b,2)  =  0.  In  addition,  no  through  thickness  distribution  of  the 
stresses  was  available.  An  approximate  approach  was  then  considered 
by  Tang  [19]  following  the  isotropic  theory  of  Reiss  and  Locke  [20]. 

The  interior  domain  (regions  removed  from  the  free  edges)  v/as  assumed 
to  be  in  a  state  of  plane  stress,  the  axial  displacement  u  was  assumed 
to  be  a  function  of  x  only,  and  the  displacement  components,  v  and  w, 
were  both  assumed  to  vanish  identically.  The  boundary  layer  equilib¬ 
rium  equations  coupled  with  the  compatibility  equations  were  split  into 
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two  problems.  Namely,  a  modified  torsion  problem  and  a  modified 
plane  strain  problem.  The  resulting  fourth  order  differential  equa¬ 
tions  were  solved  by  asymptotic  expansion  in  terms  of  the  ply  thickness 
h/2.  The  matching  of  the  boundary  layer  solution  with  the  interior 
domain  solution  was  satisfied  by  the  imposed  boundary  conditions  for 
the  two  problems.  The  results  for  a  [i45]5  graphite-epoxy  laminate 
showed  good  agreement  for  the  interior  regions  with  References  [61  and 
[7]  while  the  interlaminar  shear  stress  at  the  free  edge  was  found 
to  be  lower  than  the  predicted  singularity  of  Reference  [7].  The 
through-thickness  stress  distributions  showed  nonzero  shear  stresses 
Txz  and  on  the  free  surfaces  z  =  ±h  as  well  as  on  the  midplane 
z  =  0.  Also,  the  out-of-plane  normal  stress  vanished  on  both  the 
interfacial  planes  z  =  ±h/2  and  the  midplane  z  =  0.  This  is  unlike 
the  results  of  Reference  [7]  v/hich  indicated  maximum  values  of  on 
the  midplane  of  a  [0/90]5  laminate  and  on  the  interfacial  planes  of  a 
[±45]s  laminate.  Finally,  it  should  be  noted  that  the  approximate 
nature  of  the  formulation  in  Reference  [19]  did  not  satisfy  the  vanish¬ 
ing  stress  boundary  conditions  {±b,z)  =  0  and  cy  {?b,z)  =  0  for 
each  layer. 

1.2  THE  FINITE  DIFFERENCE  SOLUTION 

In  view  of  the  discussion  in  Section  1.1,  the  finite  difference 
solution  of  Reference  [7]  seems  to  serve  as  the  most  dependable  solu¬ 
tion  known  to  the  researcher.  This  is  due  to  the  fact  that  the  formu¬ 
lation  v/as  required  to  obtain  the  exact  elasticity  dilution  to  the 
problem.  However,  there  were  inherent  deficiencies  in  the  finite 


difference  procedures  as  pointed  out  by  Pagano  and  Pipes  [21].  To 
this  end,  numerous  tests  were  carried  out  by  this  author  to  examine 
the  "exactness”  of  the  solution  in  Reference  [7]  with  emphasis  on  its 
behavior  near  the  free  edge.  The  follovjing  observations  can  be  made. 

(1)  For  bidirectional  laminates  (0°  and  90”  plies),  all  stress 
free  boundary  conditions  were  satisfied  except  at  the  four  corners 
of  the  laminate  where  the  out-of-plane  normal  stress  did  not  vanish 
but  took  on  a  low  value.  Also,  the  sign  of  the  outer  layer  at  the 
exact  free  edge  was  found  to  be  inconsistent  v/ith  that  of  the  inner 
layer.  These  results  may  be  attributed  to  the  dissatisfaction  of  the 
equilibrium  equations  on  the  free  boundaries  as  can  be  seen  in  the  work 
of  Pipes  [12]. 

(2)  For  angle-ply  laminates  [±e]g,  neither  of  the  stress  free 
boundary  conditions,  ay(±b,2)  =  Tj(y(±b,z)  =  0,  was  satisfied  at  the 
intersection  of  the  interface  and  the  free  edge.  Both  Oy  and  t^y  at 
this  suggested  singularity  v/ere  of  an  erroneously  large  order  of 
magnitude — as  high  as  the  axial  stress  ax.  As  a  result,  the  inter¬ 
laminar  shear  stress  attained  a  maximum  finite  value  rather  than 
the  possible  infinity  predicted  by  Pipes  and  Pagano  [7].  Failure  to 
satisfy  the  vanishing  stress  boundary  conditions  at  the  four  corners, 
oz(±b,±h)  =  Tj(z{±b,±h)  =  0,  was  found  again.  Moreover,  the  sign 
reversal  of  stresses  which  was  found  for  the  bidirectional  laminates 
as  a  result  of  change  in  the  stacking  sequence,  was  not  observed  for 
the  angle-ply  laminates.  The  above  boundary  violation  may  be  due 

to  errors  inherent  in  the  solution  procedure  for  the' angle-ply 


laminates. 


It  may  be  concluded  that  despite  the  good  agreement  with  the 
results  of  References  [6].  [11]  and  [19]  in  regions  removed  from  the 
free  edges,  the  finite  difference  solution  yields  poor  results  near 
the  mathematical  singularity.  In  order  to  determine  the  proper  order 
of  magnitude  of  stress  intensity  near  the  singular  point,  a  more 
rigorous  analytical  solution  to  the  field  equations  must  be  obtained. 
Such  a  solution  was  described  by  Pagano  and  Pipes  [21]  as  a  "mathe¬ 
matical  nightmare." 

The  present  thesis  seeks  a  solution  which  predicts  accurate 
interlaminar  free  edge  stress  intensities  for  laminates.  Due  to  the 
above-mentioned  mathematical  complexities,  it  is  certainly  not  an  easy 
task.  As  described  in  the  preceding  section,  all  the  previous  investi¬ 
gations  show  a  common  result--the  plane  stress  lamination  theory  is 
recovered  near  the  central  plane  y  =  0  provided  the  laminate  is  suf¬ 
ficiently  wide  (b/h  »  1).  This  suggests  that  the  boundary  layer 
effect  is  directly  related  to  the  geometrical  ratio  b/h,  and  that  the 
stress  distribution  throughout  the  laminate  is  the  combination  of  the 
interior  region  solution  and  the  boundary  layer  solution.  The  method 
of  solution  employed  in  the  present  thesis  is  the  perturbation 
analyses  [2,  22]  developed  in  the  1940 's  to  solve  boundary  value  prob¬ 
lems  in  fluid  mechanics  and  extended  to  problems  in  solid  mechanics  in 
the  1950's.  The  isotropic  theory  of  Reiss  and  Locke  [20]  and  the 
anisotropic  theory  of  Tang  [19]  were  essentially  based  upon  such 
analyses.  The  main  differences  between  the  present  thesis  and  the 
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theory  of  Reference  [19]  are  summarized  as  follows.  (1)  The  present 
thesis  is  l)ased  upon  the  displacement  formulation  in  which  the  compati- 
bility  equ.itions  are  satisfied  automatically.  The  resulting  field 
equations  are  second  order  partial  differential  equations  in  terms  of 
the  displacement  functions.  Reference  [19]  was  based  upon  the  stress 
formulation  in  which  satisfaction  of  the  compatibility  equations 
resulted  in  fourth  order  partial  differential  equations  in  terms  of 
the  stress  functions.  (2)  For  the  interior  regions,  the  present 
thesis  dei^^rmines  the  three  dimensional  solution  to  the  reduced  govern¬ 
ing  equations  (h/b  h-  o)  while  satisfying  the  symmetry  and  antisymmetry 
conditions,  the  displacement  condition,  the  continuity  conditions 
and  the  vanishing  stress  boundary  conditions  on  the  top  and  bottom 
surfaces.  In  reference  [19]  the  displacement  components  v  and  w 
were  both  <>ssumed  to  vanish  identically  for  the  interior  regions  and 
the  axial  <lisplacement  u  was  assumed  to  be  a  linear  function  of  x 
alone  for  such  regions.'  (3)  For  the  boundary  layer  region,  the 
present  thc’sis  removes  mathematical  complexities  by  considering  the 
free  body  diagram  of  an  infinitesimally  thin  slice  containing  the 
interfacial  plane.  Such  a  limiting  analysis  provides  sufficiently 
accurate  determination  of  the  coefficients  of  the  boundary  layer 
solution  for  h/b  «  1.  The  physical  validity  of  the  composite  solution 
(interior  and  boundary  layer  solutions  combined)  is  insured  by  the 
following  requirements.  The  material  immediately  adjacent  to  the 
intcrfacial  plane  must  satisfy  the  "stretched"  governing  differential 
equations,  the  matching  principle  of  perturbation  theory  (Section  1.3), 
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the  force  and  nwment  balance  with  the  stress  resultants  on  the  central 
plane  (y  =  0),  the  self-equilibrating  conditions  of  the  out-of-plane 
normal  stress  resultant,  and  the  free  edge  stress  boundary  conditions. 
Reference  [19]  considered  two  separate  problems  for  the  boundary  layer 
region  according  to  the  even  and  odd  nature  of  the  stress  components. 

A  stress  function  following  the  isotropic  torsion  problem  was  assumed 
for  the  T°  (modified  torsion)  problem.  A  similar  function  was  then  . 
chosen  as  the  particular  solution  to  the  fourth  order  equation  of  the 
P®  (modified  plane  strain)  problem.  This  particular  solution  along 
with  the  homogeneous  solution  (5^^  order  polynomial)  constituted  the 
solution  to  this  problem.  The  combination  of  the  T°  problem  and  the 
P°  problem  failed  to  satisfy  some  of  the  stress  boundary  conditions  at 
the  free  edge  and  on  the  free  surfaces.  And  the  approximate  nature 
made  the  through-thickness  stress  distributions  incapable  of  properly 
describing  the  force  and  moment  equilibrium  and  the  self  equilibrating 
condition  at  any  level  of  z. 

1.3  PERTURBATIOfi  METHOD 

Consider  the  differential  equation 

e*y"  -  y'+y=0  ,  0ix<_l 

y(0)  =  a  ,  y(l)  =  B 

where  0  <  e  «  1 . 

Assuming  the  exact  solution  to  the  problem  is  not  available,  the 
following  approximate  steps  must  be  taken: 


(1.1) 

(1.2) 


As  t  vanishes.  Equation  (1-1)  reduces  to 

y'  -  y  -  0  (1-3) 

which  has  a  solution  of  the  form 

yo  =  ae’^  (1.4) 

where  the  superscript  "  denotes  the  solution  corresponding  to  e  =  0 
and  a  is  an  unknown  coefficient.  Solution  (1.4)  can  satisfy  only  one 
of  the  boundary  conditions  (1.2).  For  the  other  boundary  condition 
to  be  satisfied,  a  stretching  transformation  is  introduced  in  the  form 


5  =  |B  -  xl/c^  (1*5) 

where  x  >  0  and  B  is  the  boundary  limit  of  the  stretched  end  (0  or  1 
in  the  present  problem).  It  will  be  shov/n  that  this  transformation 
magnifies  a  small  region  called  the  boundary  layer  in  which  y  changes 
rapidly  in  order  to  retrieve  the  dropped  boundary  condition  at  the  end 
X  =  B.  Solution  to  the  boundary  layer  equation  must  match  the  solution 
of  the  reduced  equation  (1.3)  according  to  Prandtl's  matching 
principle  [2], 

lim  y“  ^  lim  y®*-  (1.6) 

X  B  C  -  » 

where  y®^  is  the  boundary  layer  solution. 

It  may  be  shown  [2]  for  the  present  problem  that  the  boundary 
layer  exists  near  the  end  x  =  1  and  the  value  of  >  in  Equation  (1.5) 


is  found  to  be  1.  Hence 
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(1.7) 


since  it  must  satisfy  the  first  of  Equations  (1.2).  Also,  the 
stretching  transformation  (1.5)  becomes 


(1.8) 


Equation  (1.8)  is  now  introduced  to  transform  the  original  Equation 
(1.1)  into 


^  ^  =  0 
dc2  dc 


for  e  «  1. 

Equation  (1.9)  has  the  solution 


(1.9) 


y®*-  =  c  +  de'^  (1.10) 

which  should  satisfy  the  second  of  Equation  (1.2).  Hence, 

c  +  d  =  B  (1.11) 

The  matching  principle  (1.6)  is  now  applied  as 


or. 


lim  y"  =  lim  y^l- 
X  1  ^  a. 


c  =  oe 


(1.12) 


0.^3) 


Hence  from  Equation  (1.11), 


d  =  6  -  ae 


(1.14) 
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which  yields 


y^*-  =  ae  +  (8  -  ae)e"^ 


(1.15) 


Finally  a  uniformly  valid  solution  is  fortr,ed  accordi 


ng  to  the 


equation 


yc  =  y”  +  yBL  .  (yOjBL 


(1.16) 


where  y^.  is  the  composite  solution  and  (y°)‘^‘-  represent 
part  contained  in  both  solutions. 

It  is  clear  that  in  the  present  problem 

(y“)^‘-  =  lim  yBL  =  Tin,  y“  =  , 

C  «  X  -*•  1 

hence  the  composite  solution  to  the  original  equation  i: 

yc  =  ae^  +  (b  -  ae)e'^ 


represents  the  common 


(1.17) 


(1.18) 


The  above  derivation  was  required  for  the  zeroth  order  problem 
of  Equation  (1.1).  Fcr  a  very  small  e.  the  zeroth  order  composite 
solution  (1.18)  provides  sufficient  accuracy.  For  a  relatively  larger 
e,  solution  to  higher  orders  must  be  carried  out  to  achieve  better 
accuracy.  This  is  shown  in  the  following  steps. 

The  solution  to  the  original  equation  (1.1)  may  be  expressed  as 
an  asymptotic  expansion  of  the  form 


y  =  L  c"y„(x)  ,  e  «  1 


(1.19) 


Substituting  (1.19)  into  Equation  (1.1)  results 
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J:  (c"‘^V"n  -  e"y‘n  +  =  0  (1-20) 

n=0 

Since  this  is  an  identity  equation  in  the  nonzero  parameter  e, 
the  coefficients  corresponding  to  each  n  must  vanish  for  all  x  in  the 


domain  specified  by  (1.1).  Hence, 

y’o  -  ^0  =  0  (1-21) 

y’n  -  yn  =  yn-i  (1.22) 

Also,  substituting  (1.19)  into  the  boundary  condition  (1.2)  leads  to 

yo(0)  =  «  (1.23) 

yo(l)  =  B  (1.24) 

yn(0)=ynn)=0  n>0  (1.25) 


It  is  clear  that  the  zeroth  order  problem  is  defined  by  Equa¬ 
tions  (1.21),  (1.23)  and  (1.24).  The  composite  solution  to  this 
problem  can  be  shown  to  be  identical  to  (1.18).  Also,  it  is  seen  that 
at  any  level  of  approximation  n,  yf,_|  is  known,  hence  y^  for  any  n  is 
given  by  the  first-order  equation  (1.22).  Therefore,  the  stretching 
transformation  (1.8)  should  be  continually  introduced  near  the  end 
X  =  1  where  the  boundary  condition  is  dropped.  If  the  asymptotic 
expansion 

«o 

y  =  T  e"y  (5)  ,  e  «  1 

n=0  " 


(1.26) 
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is  assumed,  the  transformed  equations  are 


«*yo 

dc2  "dc  " 


^  dyn 
dc2  '3?  Vl 


n  >  1 


And  the  boundary  condition  at  x  =  1  becomes 


(1.27) 


(1.28) 


yo(e  =  0)  =  6  (1.29) 

y^U  =  0)  =  0  n  >  1  (1.30) 

At  this  point  it  must  be  noted  that  PrandtTs  matching  principle 
(1.6)  fails  to  match  expansions  containing  higher-order  solutions. 
Instead,  Van  Dyke's  matching  principle  [2,  22]  should  be  employed  to 
obtain  a  composite  solution.  For  simplicity  in  the  present  thesis, 
only  the  zeroth  order  problem  will  be  considered,  hence  no  elaboration 
will  be  given.  Nevertheless,  it  may  well  be  an  intriguing  topic  of 
future  study. 


Chapter  II 
PROBLEM  FORMULATION 


Figure  1  shows  a  balanced  symmetric  laminate  of  2m  plies  of 
homogeneous  orthotropic  lamina  oriented  at  angles 
with  the  X  axis.  The  laminate  thickness  is  small  compared  to  other 
dimensions,  i.e.,  the  length  dimensions  are  of  the  order  L  >  b  >>  h. 
One  of  the  orthotropic  axes  of  the  laminate  coincides  with  the  z  axis. 
The  laminate  is  subjected  to  a  constant  inplane  axial  strain  As¬ 
suming  elastic  response  exists  everywhere  throughout  the  laminate,  the 
field  equations  can  be  derived  as  indicated  in  the  following  section. 


2.1  GOVERNING  FIELD  EQUATIONS 

Introducing  a  rotational  transformation  (Reference  [1])  to  the 
layerwise  orthotropic  material  leads  to  the  following  constitutive 
equations  with  respect  to  the  reference  coordinate  axes  xyz 


°x  ’ 

(k) 

‘^ll 

^12 

^13 

0 

0 

^16 

(k) 

"y 

^12 

^22 

^23 

0 

0 

C26 

^y 

^13 

^23 

^33 

0 

0 

^36 

0 

0 

0 

C44 

C45 

0 

V 

^X2 

0 

0 

0 

C45 

C55 

0 

■^xz 

Jxy. 

Cl6 

^26 

^36 

0 

0 

^66. 

-■^xy- 

where  the  superscript  k  denotes  the  k^^  layer  in  the  laminate  . 
strain-displacement  relations  in  each  layer  are 


(5^.1) 


The 
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=  u,x 
Cy  =  v,y 

Cz  =  W,Z 

(2.2) 

Tyz  -  W,y  +  v.z 
Yxz  =  W,X  +  U,2 
Yxy  =  V.X  +  u,y 

where  a  comma  denotes  partial  differentiation. 

Since  the  long  laminate  is  loaded  only  at  its  ends  x  =  tL. 

Saint  Venant's  principle  [23]  can  be  invoked  such  that  the  stresses  in 
regions  far  away  from  the  ends  are  independent  of  x.  Thus,  the 
equilibrium  equations  in  such  regions  reduce  to 


iliiZ  +  iLxz 

3y  3z 


=  0 


- 


^ 

ay  32 


=  0 


(2.3) 


Combining  equations  (2.1)  and  {2.2),  and  integrating  the 
resulting  stress-displacement  relations  (independent  of  x)  results  in 
the  following  displacement  fields  for  each  layer. 


u  =  (C^y  +  C22  +  C3)x  +  U(y,z) 

V  =  (C42  +  C5)x  -  Ci|-  +  V(y,z)  (2.4) 

w  =  -C4xy  +  Cgx  -  C2^  +  W(y,2) 

where  through  C5  are  unknown  constants  and  U,  V  and  W  are  unknown 
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functions  of  y  and  z  only. 

The  following  symmet-'y  and  antisymmetry  conditions  must  be 
imposed: 


u(x,y,z)  =  u(x,y,-z) 
v(x,y,z)  =  v(x,y,-z) 
w(x,y,z)  =  -w(x,y,-z) 
v{x,y,z)  =  -v(x,-y,z) 
w(x,y,z)  =  w(x,-y,z) 


(2.5) 


and  the  experimentally  verified  [13]  condition 

u(0,y,h)  =  -u(0,-y,h)  (2.6) 


is  imposed. 

Equation  (2.6)  leads  to  a  more  general  antisymmetry  condition 


u(o.y,z)  =  -u(o,-y.z) 


(2.7) 


for  continuity  conside.-ation. 

At  this  point,  the  even  and  odd  nature  of  the  displacements  u,  v, 
w  in  relation  to  y  and  z  can  readily  be  seen.  Substitution  of  Equa¬ 
tions  (2.4)  into  Equations  (2.5)  and  (2.7)  results  in 

Cl  =  C2  =  C4  =  C5  =  C5  =  0  (2.8) 


and 


U(y,z)  =  U(y,-z).  V(y.z)  =  V(y,-z),  W(y.z)  =  -W{y,-z) 
U(y,z)  =  -U(-y,z),  V(y,z)  =  -V(-y,z),  W(y.z)  =  W(-y,z) 


(2.9) 
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This  greatly  reduces  the  layerwise  displacement  field  functions  (2.4) 
to 

u  =  C3X  +  U{y,z) 

V  =  V{y.z)  (2.10) 

w  =  W{y,z) 

As  defined  in  the  beginning  of  the  present  chapter,  the  laminate 
is  subjected  to  a  uniform  axial  strain.  Hence  the  constant  C3  in 
Equations  (2.10)  is  nothing  but  the  applied  strain  e^. 

Combining  Equations  (2.1),  (2.2),  (2.3)  and  (2.10)  results  in 
the  following  set  of  simultaneous  partial  differential  equations  with¬ 
in  each  layer. 

C66^»yy  +  C55II.ZZ  +  C25V,yy  +  C/^^y,zz  +  (C35  +  C45)W,yz  =  0 

C26U.yy  +  C45U,zz  +  C22V,yy  +  C44V,zz  +  (C23  +  C44)W,yz  =  0  (2.11) 

(C45  +  C35)ll,yz  +  (C44  +  C23)V,yz  +  C^/^U,yy  +  C33W,zz  =  0 

The  appropriate  traction  free  boundary  conditions  are  (Fig.  1) 

Oy(tb,z)  =  0 

Txy(±b,z)  =  0  (2.12) 

Ty2(±b,z)  =0 

along  the  free  edges,  and  ' 

c^zCy.ih)  =  0 

TxzCy.ih)  =  0  (2.13) 

Tyzfy.lh)  =0 
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on  the  top  and  bottom  surfaces  of  the  laminate. 

Equations  (2.12)  and  (2.13)  may  be  expressed  in  terms  of  the 
unknown  functions  U,  V,  W  in  the  form 

{Ci2ex  +  C22V,y(±b,z)  +  C23W.z(±b,2)  +  C26U,y(±b,z)}^''^  =  0 
{Ci6tx  +  C26V,y(±b,z)  +  C35W,z(±b,z)  +  C56U,y(±b,z)  =  0  (2.12) 

{C44V,z(±b,z)  +  C44W,y(±b,z)  +  C45U,z(±bz)}^*^)  =  0 

{Ci3£x  +  C23V,y(y,±h)  +  C33W.z(y,ih)  +  C35U,y(y,ih) ^  =  0 
{C44V,z(y,±h)  +  C44W,y(y,±h)  +  C45U,z(y,±h) )  =  0  (2.13) 

{C45V,z(y,±h)  +  C45W,y(y,±h)  +  C55U,z(y,±h)}^^^  =  O 

where  the  superscripts  k  and  1  denote  the  layer  and  the  outermost 
layer  (Fig.  1),  respectively.  Equations  (2.9)  also  yield  the  following 
restrictions  on  the  displacement  fields 

{U,2(y,0)}^'"^  =  0 

{V,2(y,0))^"’^  =  0  (2.14) 

{W(y,0)}^'"^  =  0 

along  the  midplane  and 

W.z)}^*^^  =  0 

{V(O.Z)}^''^  =  0  (2.15) 

{W.y(O.z)}^''^  =  0 

along  the  central  plane  and  the  superscript  m  denotes  the  layer 
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adjacent  to  the  midplane  (Fig.  1). 

Equations  (2.11)  along  with  Equations  (2.12)  -  (2.15)  represent 
a  well  defined  boundary  value  problem.  Also,  it  is  important  to  note 
that  these  field  equations  were  derived  for  individual  layers.  Hence 
continuity  in  displacements  and  tractions  across  the  interfaces  must 
be  enforced  to  insure  completeness  of  the  solution. 

Equations  (2.11)  -  (2.15)  can  be  put  in  the  dimensionless  forms 


{Q66(h/b)2u,YY  +  QssU.ZZ  +  Q26(h/b)2v,YY  +  Q45V,ZZ 
+  (Q36  +  Q45)(h/b)W,YZ)^'^)  =  0 

{Q26(h/b)2u.YY  +  Q45U,ZZ  +  Q22(h/b)2v.YY  +  Q^^V.ZZ 
+  (Q23  +  Q44)(h/b)W,YZ)}^'^^  =  0 


^^^45  +  Q36)(h/b)U.YZ  +  (Q44  +  Q23)(h/b)V,YZ 
+  Q44(h/b)2w,yY  +  Q33W,ZZ}^'^5  =  0 

^QlZ^x  +  ^.Y(±1,Z)  +  ^,Z(±1,Z)  +  ^,Y(±1.Z)}^'^^  =  0 

^Q16^x  ^.Y(±1,Z)  +  ^,Z(±1.Z)  +  ^^').Y(±1,Z)}^''^  =  0  (2.17) 

{^.Z(±1.Z)  +  ^,Y(±1,Z)  +  ^^-U,Z(±1.Z)}‘'')  =  0 


^13^x  "  "b 


+  -Pv.Y(Y,il)  +  ^^,Z(Y,i1)  +  !^,Y(Y,‘1)’/^^  =  0 


{^.Z(Y.±1)  +  ^,Y(Y,±1)  +  ^.Z(Y.il)}^^'  =  0  (2.18) 

|^.Z(Y,il)  +  ^,Y(Y.±1)  +  ?5iy.Z(Y.*l)p^  =  0 
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{U,Z(Y,0)}^'^)  =  0 
{V,Z(Y,0)}^'"^  =  0 
{W{Y,0)}^'"^  =  0 

{U(0,Z)}^'^^  =  0 
{V(0,Z)}^'^^  =  0 
{W,Y(O.Z)}^’^^  =  0 


(2.19) 


(2.20) 


where  with  being  the  largest  stiffness  coef¬ 

ficient  of  the  ktf’  layer,  Y  =  f  and  Z  =  f  ,  the  dimensionless  co- 

D  h 

ordinates,  and  U,  V,  W  and  their  partial  derivatives  being  dimension¬ 
less  quantities. 

2.2  EQUILIBRIUM  CONSIDERATIONS 

Before  developing  the  solution  procedures,  it  will  be  shown  that 
a  close  examination  of  the  force  and  moment  equilibrium  of  a  section 
of  the  laminate  will  lead  to  significant  reduction  in  the  mathematical 
complexity.  Consider  the  free  body  diagram  ’n  Fig.  2.  Let 


tk  =  r 


k  =  1 ,2,. . . ,m 


(2.21) 


hence. 


r  t..  =  1  . 


(2.22) 


The  force  and  moment  equilibrium  per  unit  length  require 


I  Fy  =  0 


NOTE  :  Txy  ,rx7  (NOT SHOWN) 


FIGURE  2  FREE  BODY  DIAGRAM  OF  QUARTER  YZ- PLANE 
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FIGURE  3.  PARTIAL  FREE  BODY  DIAGRAM  OFQUARTER 
SECTION. 
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which  gives 

f'  <,„(0,Z)hdZ-  f'  (2-23) 

'2k 

and  Z  =  0  which  gives 

f  ay(0,Z)h^(Z  -  Z,^)dZ  =  [  02^^^(Y,Z|^)b2YdY  (2.24) 

iZk  ^0 

where  Zl  =  1  -  E  tj  is  the  elevation  of  the  interface  in  the 

j=l 

first  quadrant. 

On  the  other  hand,  the  force  equilibrium  per  unit  length 
requires 

.  r  Fx  =  0 


which  yields 


1  0 

Txy(O.Z)hdZ  =  Txz(Y.Z^)bdY 

2k  ^  'o 


(2.25) 


At  this  point,  an  important  premise  must  be  recognized  in  the 
solution  method  of  the  present  thesis.  It  has  been  numerically 
observed  in  [24]  and  examined  in  the  present  study  that  the  central 
plane  stresses  Oy^‘'^(0,Z)  and  Txy^''ho,Z)  are  essentially  constant 
in  each  layer  for  h/b  «  1  (Figs.  2,  3).  Hence  it  may  be  expressed 


mathematically  that 


(0,z)t^  =  0 


(2.26) 


E  F, 


m 

l 

k=l 


Oy 


(k) 


and 


E  F. 


m 

l 

k=l 


^xy 


=  0 


(2.27) 


for  equilibrium  considerations. 

Equations  (2.23)  through  (2.27)  together  serve  as  an  important 
tool  to  reduce  the  mathematical  complexities  in  the  present  thesis. 
Since  the  material  on  either  side  of  the  interface  Z  =  Z|^  must 
satisfy  the  governing  equations  (2.16)  and  the  boundary  conditions 
(2.17)  and  (2.20),  and  since  the  interlaminar  stress  distributions  are 
of  primary  interests,  the  boundary  layer  equations  will  be  solved  by 
considering  only  the  infinitesimally  thin  free  body  diagram  about  this 
interface  (Fig.  4).  By  doing  so,  the  boundary  value  problem  is  re¬ 
placed  by  the  free  body  force  and  moment  system  of  Fig.  4.  Thus,  the 
stress  boundary  conditions  on  the  top  and  bottom  surfaces,  (2.18),  can 
be  safely  ignored.  This  will  be  elaborated  upon  in  Subsection  2.3.2. 

2.3  PERTURBATION  SOLUTION 

As  described  in  Section  1.3,  tv;o  regions  will  be  considered 
separately.  Namely,  the  interior  region  where  the  solution  to  the 
reduced  equations  (c  0)  satisfies  boundary  conditions  at  one  end,  and 
the  boundary-layer  region  where  solution  to  the  “stretched"  equatiot'is 
satisfies  the  boundary  condition  at  the  other  end.  Matching  of  these 
two  solutions  must  be  enforced  to  insure  uniformity  of  the  resulting 


28 


composite  solution. 

2.3.1  THE  INTERIOR  REGION 

To  seek  a  straightfoiv/ard  expansion,  let 

uC')  =  E  enUn(’^)(Y,Z) 
n=0 

=  Z  e\^'^^Y,Z)  (2.28) 

n=0 

=  Z  gV^'^^(Y,Z)  e  «  1,  k  =  l,2,3,...,m 
n=0 

where  the  small  parameter  e  represents  the  geometrical  ratio  h/b. 

Substituting  these  expansions  into  Equations  (2.16)  and  equating 
coefficients  of  equal  powers  of  c  to  zero  result  in  the  following  sets 
of  equations: 

=  {<!55''o,ZZ  *  045''o,2z}‘''’  '  « 

(Q45U0.ZZ  +  Q44''o,Zzj"‘’  =  0  (2-29) 

{l)33"o.Z2r>  =  0 

e'  :  IqssU]  ,2z  +  Q45V]  _zz  ♦  (P36  *  <l45)“o.Vz}*''*  ‘  9 

{Q45>'i,ZZ  *  '!44Vi,zz  +  (Q23  +  Q44)“o,Y2}‘''’  '  “  <2-91» 

|(Q45  *  Q36)9o,VZ  *■  W44  *  923)''o,YZ  *  933W1  _zz}^'‘'  °  9 


29 


e*"  :  {Q66V2,YY  *  ^SS^r.ZZ  *  ^ZsVz.YY  *  %V,7Z 

+  (Q36  +  Q45)Wr-l.Yz}^'^^  =  0 

{Q26^r-2,YY  ^  Q45'^r,ZZ  *  Q22V2,YY  *  Q44''r,ZZ 
+  (Q23  ^  ^44^Vl,Yz}^'^^  ^ 

|(Q45  +  Q36^^r-1,YZ  *  ^^44  ■*■  ^23^'^r-l  ,YZ 

+  Q44Wr.2,YY  ■*■  Q33Wr,Zz}^’^^  =  0 

Now  the  displacement  conditions  (2.19)  and  (2.20)  give 

=  0 

=  0 

=  0  n  =  0,1 ,2, - 

=  0 

=  0 

Wn(!^)(0.Z)  =  0  n  =  0,.,2,.... 

k  =  1 ,2,. . . ,m 


and 


(m) 

Un.Z  tY.O) 

Vnjz^Y.O) 

Wn("’^(Y,0) 


Un^'^^O.Z) 

Vn^^'^O.Z) 


(2.31) 


(2.32) 


(2.33) 


Recognizing  that  the  boundary  layer  regions  exist  near  Y  =  <  1, 
the  stress  free  boundary  conditions  (2.17)  are  dropped  for  this 
interior  region. 

The  stress  boundary  conditions  on  the  top  and  bottom  surfaces. 


Equations  (2.18),  yield 
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{Qi3^x  +  P  Vn  y(Y,±1)  +  Wn  2(Y,±1)  +  ^  Un,Y(Y.il)}^^ ^  =  0 

^  Wn  j{Y.±l)  +  ^  =  0 

^  Vn,z(Y,±l)  +  Wn,Y(y.i1)  +  ^  Un.z(Y,±1)}f^)  =  O 

n  =  0,1,2 . 

For  the  lamination  theory,  a  =  t  =  t  =  n 

2  xz  yz 

(k  -  2,3,4,.,.,m)  must  also  hold  for  the  interior  region.  Hence 
Equations  (2.34)  may  be  generalized  to 

{QiS'^x  +^V„^y(Y,tl)  +%w„^z{Y>^l)  +  “b-Un.yfY.il)}^*^^  =  0 

Vn^z(Y.±l)  +  ^  Wn.Y^Y.il)  +  ^  zCY.tl =  0  (; 

+^Wn^y{Y,±l)  +  ^  =  0 

n  =  0,1,2 . 

k  =  l,2,...,m 


The  derived  symmetry  and  antisymmetry  conditions  (2.9)  lead  to 

Un^'^V.Z)  =  Unf’')(Y,-Z) 

Vn^'')(Y,Z)  =  Vn^''^{Y,-Z) 

Wn^''^(Y,Z)  =  -Wn^'')(Y,-Z) 

Unf^)(Y,Z)  =  -'J^^'^J(-Y,Z) 

Vn^'')(Y,Z)  =  -Vnf‘'){-Y,Z) 

Wnf')(Y.Z)  =  Wn(‘')(-Y,Z) 

n  =  0,1,2,.... 

k  *■  1 ,2 , .  •  a .  ,in 


(2.34) 


(2.35) 


.1. 


Equations  (2.35)  may  be  put  in  the  following  form  with  respect 
to  the  order  of  c,  similar  to  Equations  (2.29)  through  (2.31). 


(2.37) 


(2.38) 


=  ^13^''^  ^xh  +  Q33^‘'^Wo,Z^''^{y.i^)  =  0 

«45‘'’''o.z(’''f )  *  ■  « 

(k) 

••  Q23^'^^Vl.y(Y*^^)  ^  Q33^''^Wr.Z^''^y*^l) 
+  Q36^'^^U^!i!y(Y.^1)  =  0 

+  044''')w^!fjY(Y,±l)  =  0 

+  Q45^'^Vi1y'^’-1)  =  0 


Thus,  the  interior  region  problem  is  redefined  by  the  infinite 
sets  of  equations  with  respect  to  the  order  of  e. 

The  Zeroth  Order  Problem; 

Equations  (2.29),  (2.32),  (2.33),  (2.36)  and  (2.37),  (n  =  0) 
constitute  the  zeroth  order  problem  for  the  interior  region. 

The  solutions  to  Equations  (2.29)  have  the  form 


Uo^*"^  =  A^,^''^(Y)Z  +  Bq^'^'^Y) 
Vo^'^^  =  Co^‘^)(Y)Z  +  Do^'^)(Y) 
W(,(k)  =  Eof')(Y)Z  +  Fo('')(Y) 


(2.39) 


k  =  1 ,2,. . . ,m 


where  Aq^''^(Y)  through  Fq^^^(Y)  are  unknov/n  functions.  It  may  be 
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noted  that  the  form  of  Solution  (2.39)  is  similar  to  Pagano's  approxi¬ 
mate  solution  of  Reference  [17]. 

From  Equations  (2.36)  with  n  =  0,  it  may  be  shown  that 


(k) 


(k); 


(Y)  =  Cq'"'(Y)  =  Fo'"'(Y)  =  0 


(2.40) 


The  first  of  Equation  (2.37)  then  leads  to 


Eo^''^(Y)  =  - 


Qi3^'‘^Gx^ 


hence 


„  (k)  . 

° - ' 


C33 


The  last  two  of  Equations  (2.37)  are  identically  satisfied. 
From  Equations  (2.32)  with  n  =  0 

Bo^'^^O)  =  =  0  k  =  1.2,....m 


(2.41) 


(2.42) 


Also,  from  Equations  (2.36)  with  n  =  0 


Bo^''^(-Y)  =  -  Bo^'^^Y) 

Dq^'^U-Y)  =  -  Do^'^^Y)  k  =  l,2,. 


.  ,m 


(2.43) 


Equations  (2.26)  and  (2.27)  may  now  be  expressed  in  the  form 


JVxb 


^33 

(k) 
k=l 


k=l 


+  5;  C26'"'h;^Bo 


(k) 


(Y)  =  0 


(2.44) 
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^  J,  ':(i5('')hkB„'<'‘)(Y)  =  0 


(2.45) 


where  ^(Y)  and  are  the  first  derivatives  of  the  cor¬ 

responding  functions.  Note  that  the  higher  order  terms  were  neglected 
in  these  equations. 

Equation  (2.41)  implies  that  continuity  in  the  displacement 
W(Y,Z)  can  be  insured  only  when  higher  order  terms  are  included  since 
Qij  (k  =  l,2,...,m)  are  different  in  general. 

Enforcing  displacement  continuity  in  U(Y,Z)  and  V(Y,2)  results 


=  Bo^^^(Y)  =  ... 

...  =  6„<">(Y) 

(2.46) 

^^)(Y)  =  Do^2^(Y)  =  ... 

(2.47) 

Integrating  Equations  (2.44)  and  (2.45).  making  use  of  Equations 

(2.42)  and  combining  fie  resulting  Equations  with  Equations  (2.46)  and 
(2.47)  lead  to 

- >  bY 


k  =  1 ,2,... ,m  (2.48) 


D„>'''(Y)  .  .  p''5  '  . 


where 
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m  r  *^23^131 


lc=l 

m 

=  z  C 

k=1 


m  r 

fJCi2 


<^33  J 


m 


(k), 


‘’a  ■  /,''26  "k 
k=i 


^4 


- 

k=l^ 


^26^13 


16 


-33 


(k) 


rr 

%  "  ^k 


(2.49) 


As  mentioned  in  the  preceding  section,  higher-order  approxima¬ 
tions  are  not  pursued  in  this  thesis  for  simplicity.  Hence,  the 
interior  region  solutions  are  found  to  be 

jjCk)  =  +  0(c) 

yW  =  D  ^'^^(Y)  +  0(c)  (2.50) 

f  h 

w{k)  =  .  J[3_^  2  +  0(c)  k=1,2,...,ni 

^33^ 

where  Bq^*^^(Y)  and  Djj^'^^(Y)  are  given  by  Equations  (2.48),  and  0(c) 
represents  the  highest  order  term  truncated  in  the  asymptotic  expan¬ 
sion. 

2-3.2  MODIFIED  ZEROTH  ORDER  INTERIOR  REGION  SOLUTION 

Solution  (2.50)  does  not  completely  satisfy  the  vanishing  stress 
condition  (2.35)  to  the  proper  order  of  c.  This  can  be  seen 

f  k  ^ 

from  Equation  (2.38)  where  the  zeroth  order  displacements  '(Y,Z) 
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and  Vq^*'^(Y,/)  were  related  to  the  undetermined  first  order  displace¬ 
ment  w/*^)(y./).  It  has  been  described  (Section  1.3)  that  solving 
higher  order  problems  requires  more  complex  mathematical  procedures 
such  as  Van  Hike's  matching  principle.  Also  in  Section  2.2,  it 
was  shown  th.ii  the  uniform  stress  distribution  on  the  central  xz  - 
plane  (Y  =  0).  a  numerical  result,  is  utilized  as  an  important  tool 
to  reduce  mait>.‘matical  complexity  for  the  boundary  layer  region. 
Therefore,  a»  improved  zeroth  order  interior  region  solution  to  evalu¬ 
ate  better  si-vss  intensity  near  the  central  plane  is  certainly  quite 
demanding. 

To  seek  Mich  an  improvement.  Equation  (2.39)  along  with  Equation 
(2.40)  are  nc»  required  to  satisfy  the  stress  conditions  (2.35) 

exactly.  Equ.rions  (2.32).  (2.33)  and  (2.36)  remain  satisfied.  The 
following  equ.-.ions  are  obtained. 


W, 


=  r  (k) 


0  =  ^0 
=  p  (k) 


V  f')  =  D 
0  *^0 


(k) 


Z 

(Y) 

(Y) 


‘^13S  b 


^23 


Do*(Y)  +  ^Eq  +  ^Bo'(Y)} 


(k) 


Hhere  Eq^*^)  nov  -ecomes  an  unknown  constant  for  the  layer. 
Again  er  •  cing  contiiuity  in  displacements  Uq^*^^(Y,Z)  and 
rc- •'Actively,  yields 


(2.51) 

(2.52) 

(2.53) 

(2.54) 


(Y)  =  iro(Y) 


(2.55) 
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and  =  Do^^^(Y)  =  ....  =  Do^'')(Y)  =  Do(Y)  (2.56) 

Hence,  Equation  (2.54)  may  be  written  as 

r  W  r  r 

? - Do'(Y)  +  - Bo*(Y)  =  0  (2.57) 

(k) 

Continuity  in  the  displacement  (Y,Z),  as  developed  in  Sub¬ 
section  2.3.1  (Equation  (2.41)),  will  be  insured  only  by  higher  order 
considerations,  hence  it  is  not  imposed  as  a  physical  requirement  in 
the  present  modification. 

Finally,  recalling  Equations  (2.26)  and  (2.27)  gives 


"  m 

L  ,  ' 

(k)  ' 

Z 

k=l 

C12  IT  ^0 

i 

•’k 

* 

(k) 


^V-^Do‘(Y) 


k=1 


j 


and 


Bo‘(Y)  =  0 


m 

l 

k=l 


L  ^36  ^ 

(k)  1 

f 

C16  +  IT  ^0 

.  J 

f’k 

H  * 

n.  r 

m  C2-  h 


k=l 


-jffo'W 


(2.58) 


k=i  " 


B„’(Y)  =  0 


(2.59) 


Since  there  are  m  +  2  simultaneous  equations  ((2,57)  -  (2.59)) 
for  the  m  +  2  unknowns  Bo'(Y),  Dq’(Y),  and  Eq^*^\  the  modified  zeroth 
order  interior  solution  can  be  readily  determined. 

To  show  the  improvement  made  in  the  present  modified  Interior 
region  solution,  two  numerical  examples  are  given  in  Tables  1  and  2. 


TABLE  1* 


[0/90]s  (c  =  0.133) 

Displacement 

Finite  Difference  ZIR**  Solution 

Modified 
ZIR  Solution 

uf^)/{exbY) 

0  0 

0 

v('')/(exbY) 

-0.0397  -0.0391 

-0.0396 

wt^V(cj^hZ) 

-0.2467  -0.2534 

-0.2448 

w(2)/(exhZ) 

-0.2055  -0.2172 

-0.2072 

*  Material:  graphite-epoxy  laminate  with  constant  ply 
thickness. 

**  Zero^^ 

order  interior  region. 

TABLE 

[45/-45]s  (c 

2 

=  0.133) 

Modified 

Displacement 

Finite  Difference 

ZIR  Solution 

ZIR  Solution 

U^'^)/(exbY) 

0 

0 

0 

-0.7409 

-0.7298 

-0.7433 

Mn)/(e^hZ) 

-0.0607 

-0.2354 

-0.0604 

V(cxbZ) 

-0.0613 

-0.2354 

-0.0604 

As  expected,  the  Modified  ZIR  solution  yields  jnore  reliable 
results  than  the  ZIR  solution. 
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Hence  the  former  will  be  employed  throughout  the  present  thesis 
to  evaluate  the  central  plane  stress  intensity  for  the  determination 
of  the  unknown  coefficients  in  the  boundary  layer  region  solution. 

2.3.3  THE  BOUNDARY  LAYER  REGION 

Consider  the  first  quadrant  of  the  yz  -  plane  as  shown  in 
Figure  5.  Introducing  the  stretching  transformation 

n  .  (2.60) 

near  the  free  edge  Y  =  1  to  the  governing  equations  (2.16)  results  in 
the  following  equations  for  this  quarter  plane  of  the  laminate. 

Ikl 

^%6^*nn  ^  ^55*^*ZZ  *  *  ^'^’ZZ  "  ^^36  ^  Q455''^*riZJ  "  0 

(k) 

^^26'^*nn  ^AS'^’ZZ  ^  *  Q44V»ZZ  '  ^^23  ^  044)'^’nZ^  ""  ° 

(1^)  . 

t-(Q45  +  Q36)U»nZ  ‘  ^^44  Q23)''’nZ  Q44'''’nn  ^33W>ZZ>  "  ^ 

(2.61) 

To  satisfy  Prandtl's  matching  principle  (Section  1.3),  assume  the 
following  composite  expansions 

=  [Bo(Y)  +  pQe^o'’  cos  ao  +  0(..) 

=  [Do(Y)  +  py®"'  cos  ao  Z]^*")  +  0(0  k  =  1.2,. ..,m  (2.62) 

w(k)  =  [EpZ  +  Soe^®"*  sin  a-o  +  0(e) 

where  Bo^*^^(Y).  Dq^'^^Y)  and  Eq^*^^  are  the  Modified  ZIR  solution 

fkl  (kl  (k) 

given  by  Equations  (2.57)  -  (2.59),  Pq'  O  Rq  and  Sq'  are  un- 

(k)  .  .  . 

determined  coefficients,  and  aQ  are  undetermined  positive  quantities 
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in  radians.  The  subscript  o  implies  the  zeroth  order  solution  for  the 
boundary  layer  region. 

Substituting  Equations  (2.62)  into  Equations  (2.61)  and  neglect¬ 
ing  the  0(c)  terms  results  in  the  following  set  of  three  simultaneous 
algebraic  equations  corresponding  to  the  c°  order: 

o  0  2  ( k ) 

^^^66^0  ■  %5“o  ^^0  *  ^^26^0  ‘  ^45“o  ^^o  ”  ^^36  '*■  ^45^^o“o^o^  "  ® 

^^^26^0^  "  ^45“o^^^o  ^^22^0^  "  ^44“o^^^o  "  ^^23  ^44^^o“o^o^^  ^  ® 

{(Q45  +  Q36)^o“o^o  ^^44  *^23^^o“o^o  ^  ^^44^0  '  ^33“o^^^o^  "  ® 

k  =  1,2.... .m  (2.63) 


For  each  nontrivial  term  of  Solution  (2.62)  to  exist  the 
determinants  of  these  algebraic  equations  must  vanish  individually. 
Thus , 


^66^0^ 

-  Q55“o' 

^26^0^ 

■  ^45'"o^ 

-(Q36  ^45)^0'^ 

(k) 

^26^0^ 

-  Q45“o^ 

^22^0^ 

■  ^44^0^ 

-(O23  +  Q44^^o"o 

=  0  (2.64) 

(Q45  + 

Q36)^o“o 

(Q44  + 

Q23)^o'^o 

^44^0^  "  ^33'=‘o^ 

k  =  1 ,2,...  ,m 

These  sixth  order  equations  may  be  regarded  as  third-order  equa¬ 
tions  by  the  classical  treatment  [25]  and  the  method  of  complex 
variables  [26]. 

The  six  roots  are  found  to  be  in  the  form 

{Xo(l,2)  =  ±  ? 

{Xo(3,4)  =  1  b  (2.65) 

{Ao(5,6)  =  i  C 
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— /ij>  — (k)  — (k) 

where  a  b  ,  c  '  are  three  positive  constants  in  terms  of 
material  constants  of  the  k^^  layer  (see  Section  3.2).  For  matching 
consideration,  however,  the  positive  roots  must  be  dropped  since  they 
lead  to  exponential  growths  of  the  displacements  for  large  n  (or  small 

Y). 

Upon  determining  the  characteristic  roots  from  Equations  (2.65), 
Solution  (2.62)  takes  the  following  general  form. 

uM  =  {Bo(y)  +  (P^e-^“o^  +  +  P3e^'^on)  +  0(c) 

=  {Do(Y)  +  (R^e’"“on  +  +  R3;^“onj  ,,3  2j(k)  ^ 

wf')  =  {EqZ  +  (S^e^^on  +  526^“°^  +  S3e''^^o'’)  sin  +  0(c)  (2.66) 

where  are  replaced  by  etc. 

With  the  above  solution,  the  stress  boundary  conditions  (2.17a), 
(2.17b)  and  (2.17c)  transform  to 


([Q26(3Pi  +  bP^  +  CP3)  +  Q22(aRi  +  bP2  +  CR3) 

+  Q23(Si  +  $2  +  S3)]ao  cos(fto  A) 

=  -  C(Qi2  +^Eo)c^  +^Do'(±1)  Bo'(-tl)]h)^''^ 

(2.67) 


t[Q66(aPl  +  bP2  ^  <^^*3)  Q26(«'^1  +  ‘^'^2 

+  036(^1  +  S2  +  S3)]ao  cos  (ao  Z^) 

=  -  [(Q16  +  ^fo)--x  +  -F  Do‘(^1)  +  ^Bo'(^l)h}^''^ 


(2.68) 
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{Q44[(Ri  +  R2  +  R3)  -  (S^a  +  526"  +  S3C)]  +  Q45(Pi  +  P2  +  P3)}^*'^  =  0 

k  =  l,2,...,m  (2.69) 

Note  that  the  right  hand  sides  of  Equations  (2.67)  and  (2.68)  are  all 
known  quantities  from  the  interior  problem.  Since  there  are  ten  un¬ 
known  coefficients  in  the  k^^  layer,  solving  simultaneously  three 
equations  from  the  boundary  conditions  (2.67)  through  (2.69),  and  six 
equations  from  Equations  (2.63)  leads  to  the  determination  of  the  nine 
unknown  coefficients  in  terms  of  qq.  The  accuracy  of  the  coefficients 
thus  obtained  can  be  readily  checked  by  the  self-equilibrating  condi¬ 
tion  of  the  stress  resultant, 

^  ^  j  =  0  (2.70) 

for  any  level  of  Z  (Fig.  2). 

Finally,  imposing  the  moment  equilibrium  conditions  (2.24)  and 
the  force  equilibrium  conditions  (2.23)  and  (2.25)  determines  the 
values  of  ^  and  tan  (oto^^^^Zj.)  to  their  orders  of  accuracy. 

In  summary,  the  zeroth  order  interior  solution  (ZIR)  was  ob¬ 
tained  by  letting  h/b  go  to  zero.  The  Modified  ZIR  solution  improved 
the  ZIR  solution  by  satisfying  the  vanishing  stress  boundary  condi¬ 
tions  (2.35)  exactly.  The  zeroth  order  boundary  layer  solution  was 
obtained  by  transforming  the  governing  equations  and  the  boundary  con¬ 
ditions  (2.17)  at  the  free  edge.  The  matching  principle  was  satisfied 
by  the  composite  solution,  and  the  self-equilibrating  condition  of  the 
interlaminar  normal  stress  resultant  was  employed  to  check  the 
accuracy  of  the  calculated  coefficients.  The  continuity  conditions 
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in  both  displacements  and  tractions  were  imposed.  And  the  force  and 
moment  equilibrium  of  the  composite  solution  with  the  central  plane 
stress  resultants  were  satisfied  for  the  k^h  interfacial  plane 
Z  =  Zk  (Figs.  3.  4). 


Chapter  III 
SPECIAL  LAMINATES 


The  solution  method  developed  in  the  preceding  chapter  applies 
to  balanced,  symmetric  laminates  with  variable-thickness  plies.  For 
certain  special  cases  the  field  equations  are  greatly  simplified  by 
the  vanishing  of  some  elements  in  the  stiffness  matrix.  Among  the 
various  laminates  studied  in  the  literature  (Section  1.1)  are  the  bi¬ 
directional  laminates  [0/90]5  and  [90/0]5,  and  the  angle-ply  laminates 
[0/-0]s  C-e/0]s'  These  two  laminates  will  be  considered  in  this 

chapter. 


3.1  BIDIRECTIONAL  LAMINATES  WITH  CONSTANT  PLY  THICKNESS 

When  the  orientation  of  the  fibrous  layer  is  either  O'"  or  90“ 
with  respect  to  the  x  axis,  the  constitutive  equation  reduces  to 


Ox  ' 

(k) 

cn 

Cl  2 

Cl3 

0 

0 

0  ■ 

(k) 

‘'y 

Cl2 

C22 

C23 

0 

0 

0 

"y 

OZ 

^13 

C23 

C33 

0 

0 

0 

^z 

V 

0 

0 

0 

C44 

0 

0 

^yz 

^xz 

0 

0 

0 

0 

0 

YxZ 

-■^xy. 

.  0 

0 

0 

0 

0 

C66. 

-^xy. 

k  =  1 ,2,. . . ,m. 


Consider  the  laminate  consisting  of  2m  layers  with  the  stacking 
sequence  [0/90/0/90/0/90. . .3^  as  shown  in  Fig.  6a. 


I 


I 
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Equations  (2.16)  reduce  to 

2 

-%6  [Sj  ^’YY  *  Q55  =  0 

{Q22  [if  V.yy  +  Q44  V.zz  +  (Q23  +  Q44)[^j  =  0  (3.2) 

r  ^  2 

{(Q44  +  Q23Hb]  '^’YZ  *  ^44  [b]  W,yy  +  Q33  W.zz)^*'^  =  0 

Note  that  the  first  equation  is  an  independent  equation  whose  complete 
solution  may  be  assumed  in  the  form 


"  AnY  (k) 

r  a^^  e  cos  Op  Z) 
n=0 


I  K  I 

where  a^  (n  =  0,1,2, - )  are  unknown  coefficients. 

Substituting  Equation  (3.3)  into  the  first  of  Equations  (3.2) 


yield 


»  (k)  _r  Mb  \(k) 
"(1.2)  ■  r/Oee  "  V 


n  =  0,1 ,2 . 


Hence 


|(k)  = 


Solution  (3.5)  automatically  satisfies  the  first  of  the  displace¬ 
ment  symmetry  conditions  (2.19)  and  (2.20).  The  second  equation  of  the 
free  edge  stress  boundary  conditions  (2.17)  reduces  to 


|^U,Y  (n,Z)|^'^)  -  0 


(3.6) 
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Substituting  Solution  (3.5)  into  Equation  (3.6)  results  in 


n  =  0,1 ,2 . 


(3.7) 


hence , 


(3.8) 


everywhere  in  this  laminate. 

(k)  (k) 

This  leads  to  the  vanishing  of  the  shear  stresses  and 

throughout  the  bidirectional  laminate  as  may  be  physically  expected. 

fk)  fk) 

The  modified  ZIR  solution  (Subsection  2.3.2)  for  >  and  W' 
may  be  determined  by  solving  Equations  (2.57)  -  (2.59)  simultaneously 
with  Bo'(Y)  vanishing  identically.  For  the  simplest  case  of  the  four 
layer  symmetric  [0/90]5  laminate,  the  displacements  are  found  to  be 

v(k)  =  D^(k)(Y) 

y(l)  ,,  (3.9) 

w(2)  =  e^{2)2 


where 


(2)  (')  ,  ,  f  (2)  r(')r(2)l 

c  (2)  -  ^23i33  c  '0  J.  E13  .  '>3^23  L  h 


^23  43 


^  cTriFr2T 

43  23  33 


(3.10) 


(3.11) 


.43  " 


^33^^  r  .  ^i3^*^x 


Tnr 


(3.12) 


Introducing  the  stretching  transformation 
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n  -  Y) 

'  [Sj 

into  the  remaining  coupled  equations  of  (3.2)  results  in  the  boundary 
layer  equation  in  the  form 


^^22  Q44  V.2Z  -  (Q23  +  Q44)  =  0 

{-(Q44  +  Q23)  V,^2  +  Q44  W,r,T  +  Q33  W,22}^'^^  =  0 


(3.14) 


Following  Subsection  2.3.3,  the  displacements  are  assumed  to  be 
v(*')  =  Do('^)(Y)  +  [R^e  °  cos  oq  +  0(t. ) 

(3.15 

u(^)  r  f  (1^)7  4.  r<-  _^'0^  .  ..ifH 

‘'0  Z  +  [So  e  sin  +  0(c) 

where  Dq  (Y),  Eq  ^  are  known  quantities  from  the  modified  ZIR  solu¬ 
tion.  For  the  four  ply  [0/90],  laminate,  they  are  given  by  Equations 
(3.10)  -  (3.12).  and  are  unknov/n  coefficients. 

Substituting  Equations  (3.15)  into  Equations  (3.14)  results  in 

the  algebraic  equations  for  the  zeroth  order  boundary  layer  problem  as 
follows 


{(Q22  ^0^  -  Q44  ''o^)'^o  -  (^23  +  “0  =  0 

UQ44  +  Q23)>o  ^o  Rq  ^  (Q44  ^0^  -  P33  ° 

For  a  nontrivial  solution,  the  determinant  of  these  equations 
must  vanish.  Thus. 


^22  ^0^  ■  ^44  '*0^ 

(Q44  +  023^^0  ^‘0 


-(Q23  +  Q44)>o  ''0 

^44  '0^  "  O33  ''0^ 


(3.17) 
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whence 


(k)  =  . 


.2  -  4  2s 


'22 


T7? 


(k) 


(k) 


(3.18} 


.(k)  = 


Q33  _  Q23Q23  ^^23^44] 

Q44  '  O22Q44  j 


(k) 


(3.19) 


Note  that  the  2  x  2  determinant  is  only  a  principal  minor  of  the 
determinant  in  Equations  (2.64).  For  conventional  composites,  such  as 
graphite-epoxy  and  boron-epoxy  laminates, 

(3.20) 


p2-4 


0331’''’ 


O22 


>  0 


Hence  Equation  (3.18)  yields  two  pairs  of  real  roots.  For  match¬ 
ing  considerations,  the  positive  roots  are  dropped,  and  the  zeroth 
order  composite  solution  (Section  1.3)  takes  the  following  form 


V  (k)  =  Do<‘‘>(V)  ♦  ((R,  s'®'  ♦  "2  “0  2! 


(k) 


w  +  {(S]  e^l  +  S2  “°'’)sin  iq 

-(k) 


(k) 


(3.21) 


where 


(k)  = 


,(k) 


1 - 

p  + 

|p2  -  4  2S1 
1  ‘^22] 

M2 

_ 

2 

! 

f 

P  - 

2  .  ^33 

P  -  4  - 

^22. 

1/2 

(3.22) 


(k) 


62'  '  7  2 

and  Prandtl's  matching  principle  (Section  1.3)  is  satisfied. 
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The  first  and  third  of  boundary  conditions  (2.17)  lead  to 

S2R2)  Q23(^1^2]'^o  ^02  f“o  ^k) 

=  -  cxh(Qi2  +  Q22  Do‘(Y)  +  Q23Eo)>^^^  ^3  23) 

{Q44  (Ri  +  R2  -  Si6i  -  S2B2)>^''^  =  0 

From  Equations  (3.16),  additional  relations  between  and 
R2  and  $2,  are  obtained  as  follows 

[(022^1^  -  044)^1  -  (Q23  Q44)eiSi]^''^  =  0  (3-24) 

[(02262^  -  Q44)R2  -  (Q23  +  Q44)S2S2]^''^  =  0  (3.25) 

Solving  Equations  (3.23)  -  (3.25)  simultaneously  results  in 

Rl^*^^,  R2^*'^>  $1^^^  and  $2^^^  in  terms,  of  cos  t  4)) 

where  0  <  c  «<  1  and  Zj^  is  defined  in  Equation  (2.24)  and  (2.25). 

Equations  (2.23)  and  (2.24)  then  lead  to  the  determination  of 

and  tan  Z|(  -  4)  to  their  orders  of  accuracy. 

Thus  the  complete  solution  for  the  zeroth  order  displacement 

function  are  obtained  and  the  interlaminar  stresses 

between  the  k^^  layer  and  the  (k+1)^^  layer,  Z|^  =  1  -  r  t^  ,  may  be 

r=1 

readily  calculated  from  the  strain-displacement  equations  (2.2)  and  the 
constitutive  equations  (3.1). 

3.1.1  [0/90]s  GRAPHITE-EPOXY  LAMINATE 

As  a  numerical  example,  the  four-ply  [O/QOlj  graphite-epoxy  lami¬ 
nate  with  constant  ply  thicknesses  (Fig.  7a)  is  considered.  The  stiff¬ 
ness  coefficients  (after  transformation)  are  listed  below. 


<  o 
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t _  i 

0  (1) 

^1/ 

1  A 

90  (1) 

1 

90  (2) 

1  A 

0 

^  0  (2) 

0  (I) 


AE 


(b)  QUARTER  Y-Z  PLANE 
OF  TOP  LAYER 

i  i 


o<C«<i 

(c)  FBD  OFZ=l/2C 


(d)  [90/ ol  LAMINATE 


90(1) 


T 


(e)  QUARTER  Y-Z  PLANE 
OF  TOP  LAYER 

rU 


TfT 


o<C<«i 

(?)FBD  0FZ  =  I/2C 


FIGURE  7.  FOUR  PLY  BIDIRECTIONAL  LAMINATES 
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0  (x  10"^  psi) 

c{]^  =  20.2 


0.56 

2.21 

0.56 

0.48 


cH)  = 

‘'12 

rd)  = 

22 

cip  = 

c(')  . 

S3 

ctp  .  2.3, 


rd)  = 
‘'44 

rd)  = 
55 

rd)  = 

‘'66 


0.85 

0.85 

0.85 


90  (x  10'^  psi) 

2.21 


rd) 

hi 


cd)  =  0.56 

cd)  =  20.2 

Cj|)  =  0.48 


cd)  = 

23 

cd)  = 
33 

cd)  = 

‘'44 

cd)  = 

r(2)  = 
''66 


0.56 

2.21 

O.PS 

0.85 

0.85 


(3.26) 


From  Equation  (3.8),  the  axial  displacement  function  U  vanishes 
everywhere  in  the  laminate. 

The  interior  region  solutions  (3.9)  are  found  to  be 


(2) 

=  Vq‘  '  =  -0.0396  c^bY 

(3.27) 

d)  =  -0,2448 

(3.28) 

d)  =  -0.2072  c^hZ 

(3.29) 

From  Equations  (2.2),  (3.1),  and  (3.27)  -  (3.29),  the  central 
plane  (Y  =  0)  stresses  are  found  to  be 


ayd)(o,Z)  =  0.3552  cx(10®  psi) 


(3.30) 


/ 
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=  0.8385 

=  -0.0028  «2 

R2(^^  =  -1.1776  (i-i 

=  0.0984  ^2 

=  -1.0619 

=  0.0625  ^2 

$2^^^  =  0.9298 

$2^^^  =  -0.0134  'P2 

where 

cxh 

"  ^’'7^cos(ao(^)(i 

+  0) 

A  _  X 

2  ajj(2)cos(ao(2)(^ 

-  c)) 

0  < 

c  «<  1 

The  self-equilibrating  condition  f: 

=  0,  Equation  (2.70 

be  'written  in  the  form 

-f  a2^'^^(n,  ^  ±  c)hdn  =  0 

0  <  c  <«  1 

(3.34) 


(3.35) 


(3.36) 


(3.37) 

Substituting  the  coefficients  of  Equations  (3.34)  into  Equations 
(3.21)  and  the  constitutive  equation  (3.1)  determines  the  stresses  on 
either  side  of  the  interface  Z  =  ^.  It  may  be  shov/n  that  Equation 
(3.37)  is  satisfied  identically.  This  further  confirms  the  correct¬ 
ness  of  the  calculated  coefficients  of  Equations  (3.34). 

Equations  (2.23)  and  (2.24)  now  become 

I  -  c)hdn  =  -  ht  ^(10^)  (3.38) 


i: 


-  T 


yz 


(2) 


0.3552 


(n,  j  -  Ohdn  =  -  -  f-  hc^(10") 


(3.39) 
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f  -  02(l)(n,  j  +  c)bh(1  -  en)dn  =  -  be^(10®)J-  (3.40) 

1°  -  I  -  0bh(l  -  en)dn  =  br.x{10®)j  (3.41) 

oo 

0  <  c  '«  1 


where  ^  is  the  approximate  distance  of  the  resultants 
from  the  interfacia''  plane. 

To  compare  with  the  numerical  results  of  Pipes  and  Paqano  [7], 
the  interlaminar  stresses  are  calculated  based  on  the  90°  ply  (the 


lower  layer).  Equation  (3.39)  leads  to 
tan 


'=2^  -  . 


w 


=  0.5 


0  <  c  1 


whence  a, 


(2) 


0.180,  8.9868,  15.4505, 


Equation  (3.41)  gives 


(2)  = 


2.8284 


(3.42) 

(3.43) 

(3.44) 


From  Equations  (3.38)  -  (3.41),  it  is  clear  that  the  stress  resultant 

is  of  order  0(hej(10^)  v.hile  the  couple  moment  is  of  order  0(h ). 

Hence,  requiring  exact  satisfaction  of  Equation  (3.39)  and  approximate 

(2) 

satisfaction  of  Equation  (3.41)  fixes  the  value  of  at 

=  8.9868  (3.45) 

While  Tjjy  and  vanish  throughout  the  laminate,  the  other 
stress  components  are  obtained  in  the  following  zeroth  order  fonn 


(2) 


c)  =  [2.08837  +  (0.0275 

+  0.0051  e^^‘^^^^)]f;jj(10^  psi)- 


(3.46) 
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14.158r 


d  =  -  (0.3552  +  0.0546 

-  0.4098  psi) 

Ty/2^(n,|--  c)  =  0.3865(i’-®®^'’  -  psi) 


(3.47) 


(3.48) 


/^^(n  1-  c)  =  (0.1356  -  o. 


0185  e^'®®^")cx(10®  psi)  (3.49) 


where  0  <  ^  1 

The  last  two  components  of  stress,  the  interlaminar  stresses, 
are  plotted  and  compared  with  numerical  results  in  Figs.  8  -  11. 

If  the  stacking  sequence  of  the  laminate  is  reversed  to  [90/0]5, 
(^^9-  7d) ,  the  derivation  of  these  interlaminar  stresses  is  as  indi¬ 
cated  in  the  follov/ing  subsection. 

3.1.2  [90/0]s  graphite-epoxy  LAMINATE 

While  vanishes  everywhere  in  the  laminates,  the  modified  ZIR 
solution  for  and  can  be  obtained  by  interchanging  the 
superscripts  in  Equations  (3.27)  through  (3.41). 

To  compare  with  the  numerical  results  L7]»  the  0°  ply  is  now 
used  as  the  reference  layer  for  the  interlaminar  stresses. 

The  stress  components  in  zeroth  order  forms  are  obtained  as 

-  c)  =  [20.04  -  C.074(i^-566n  +  (10®  psi)  (3.50) 


<7y^^^(n,|--  c)  =  (0.3552  +  1.5452 

_  1  onrtA  r8.104r|.  /,_6  _ 


-  1.9004  e  •  psi)  (3.51) 

1  -  d  .  7,6996(5’-5“’  .  i®' p,i)  (3,52, 

hn.2  ■  O  =  (-1.9004  1.5452  e^''^^'')£jj(10®  psi)  (3.53) 

0  <  c  1 

where  the  last  two  are  the  interlaminar  stresses. 


(!Sd)g_0l  X 


FIGURE  8.  INTERLAMINAR  SHEAR  STRESS  Tyz  AS 
A  FUNCTION  OF  €  FOR  [0/90],, 


►Y 


58 


FIGURE  9.  INTERLAMINAR  NORMAL  STRESS  O2  AS  A  FUNCTION  OF€ 
FOR  [0/90]c 


0.530 


—  —  O  Q  Q  Q 

(!Sd)gOI 
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(ISdjg-OI  X 


FIGURE  13.  INTERLAMINAR  NORMAL  STRESS  Oj  AS  A 
FUNCTION  0F€  for  [SO/Olg. 


RESULTS  FOR  Ty?  IN  [90/0] 


FOR  O',  IN  [90/0] 
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These  stresses  are  plotted  in  Figs.  12  -  15. 


3.2  ANGLE  PLY  LAMINATES  WITH  CONSTANT  PLY  THICKNESS 

For  contemporary  fiber-reinforced  composites  having  three  mutual 
ly  perpendicular  planes  of  elastic  symmetry,  vanishes.  If  the 

laminate  consists  of  one  material  with  symmetric  [eZ-ejj  or  [-o/o]^ 
orientations  (Fig.  6b),  it  is  called  an  angle  ply  laminate  and  the 
following  relations  between  material  constants  are  found, 

Cij^^  =  ,  i  =  1,2,3  and  j  =  1,2,3 

CkP^  =  CkP^  ,  k  =  4,5,6  (3.54) 

r  -  r  (2)  n  -  1  ?  * 

Cn6  ■  "  Cn6  »  n  -  1,2,3 


The  modified  ZIR  solution  gives 


Uo(l)  =  Uq(2)  =0 

V  (1)  =  V  (2)  =  .  ^^12^33  -  ^13^23)  Y 

0  0  (^22^33  "  ^3^23^^^ 


(2)  _  (Ct3C22  - 


(C 


22^33 


-  C23C23) 


X 

IT)- 


(3.55) 


On  the  central  plane  (Y  =  0),  the  stresses  are  obtained  from 
Equations  (3.55),  (2.2),  and  (3.1)  as 

ay^l)(0,Z)  =  -  Oy(2)(o,z)  =  0 


(3.56) 
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Txy^^^O.Z)  =  -  T^y^2)(o'2) 

=  fr,.  -  ^_26^^12^33  -  -  C,,C,,h(1) 

L  (C22C33  -  C23C2I) - J  ‘^x 

(3.57) 

The  first  equation  indicates  that  the  zeroth  order  solution 
(3.55)  contributes  no  transverse  normal  stress  throughout  the  angle -ply 
laminate.  For  the  laminate  to  be  in  equilibrium,  two  self-equilibrat¬ 
ing  conditions  in  addition  to  Equation  (2.70)  should  be  expected  to 
hold  (Fig.  16).  Recalling  Equations  (3.38)  through  (3.41),  the  follow¬ 
ing  equations  may  be  established. 


*  c)hdn  =  0  k  =  1 ,2 


0  <  c  <«  1 


‘L  =  0  k  =  1,2 


(3.58) 


(3.59) 


Furthermore,  the  shear  stress  resultants  and  must 

also  be  in  equilibrium  (Fig.  3  with  m  =  2)  as  indicated  in  the  follow¬ 
ing  equations. 


'L  ^xz^^^(^»2^hdn  +  R^y^^)  =  0  where  R^y^^^  =  'xy^^^(0,Z)hdZ  (3.60) 

2 

~J„  ^xz  +  Rxy^^^  =  0  where  R^^y^^)  „  T^y^^^(0,Z)hdZ  (3.61) 

The  characteristic  equation  (2.64)  leads  to  two  identical  sixth 
order  algebraic  equations  for  both  layers.  Three  positive  roots  to 
this  equation  must  be  dropped  for  matching  considerations.  The  compos¬ 
ite  solution  will  be  in  the  form  of  Equation  (2.66)  with  0(c) 
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truncated.  A  numerical  example  is  presented  in  the  following  sub¬ 
section. 

3.2.1  [45/-45]5  GRAPHITE-EPOXY  LAMINATE 

Consider  the  [45/-45]5  graphite-epoxy  laminate  of  constant  ply 
thickness  h/2  (Fig.  17a).  The  stiffness  coefficients  (after  transfor¬ 
mation)  are 


45 (x  10'^  psi) 

=  6.745 

=  5.045 
=  0.521 
C22^^  =  6.745 
€23’^  =  0.521 
Cj'"  =  2.213 

(1 )  _  c  =  -4  506 

he  26 

63^^^  =  -0.04387 

f  -  r  -  0  86 
M4  ■  %5  ■ 

CjJ"  .  5.33 
C  =  0 


-45 (x  10'6  psi) 

=  6.745 

=  5.045 
=  0.521 
C22^^  =  6.745 
0^3^^  =  0.521 
Cjf  >  =  2.213 

63^^^  =  0.04387 

r  (2)  .  r  (2)  .  .  r.; 

^4  ■  ^5  ■ 

(2)  . 

‘'66 

=  0 
45 


The  modified  ZIR  solution  (3.55)  gives 


U, 


(1)  =  u  = 


=  0 


V  (T)  =  V  =  - 
*0  *^0 


0.7433  tj^bY 


W, 


(1) 


W, 


(2)  _ 


=  -0.0604  f.yhZ 


(3.62) 


<  X 
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(b)  QUARTER  (f)  QUARTER  YZ- PLANE 

Y2-  PLANE  OF  TOP  LAYER  OF  TOP  LAYER 


0<5«<l 

(d)FBD  OFZ  =  )^2±C  (h)FBDOFZ=^2±^ 


FIGURE  17.  FOUR  LAYER  ANGLE  -  PLY  LAMINATE 


Equation  (3.57)  yields  =  -Tj^y^^(0,Z)  =  1.l54cj^(10®  psi) 

The  bcundary  layer  equations  (2.61)  through  (2.64)  yield  the 
algebraic  equation  for  both  layers 


Xq®  -  2.5460  +  1.6337  -  0.1202  =  0  (3.63) 


which  is  readily  transformed  to 

-  2.5460  J  +  1.6337  oj  -  0.1202  =  0  (3.64) 

1/2 

by  letting  Xq  -  ±  (3.65) 

Furthermore,  setting  u  =  y  -  l/3(-2.5460)  =  y  +  0.8487  (3.66) 

results  in  y^  -  0.5269  y  +  0.0438  =  0  (3.67) 

bet  y  =  Po  +  %  (3.68) 


and  substitute  it  into  Equation  (3.67).  The  resulting  set  of  algebraic 
equations  are 


Po^  +  =  -0.0438 

PoPo  =  0-^756 


(3.69) 


which  give 

=  (0.0736){cos(2k7r  +  107.3“’)  +  i  sin(2kii  +  107.3“)} 
qjj^  =  (0.0736){cos(2k:T  +  252.7“)  +  i  sin(2kTT  +  252.7“)} 

k  =  0,1,2 


By  applying  DeMoivre's  formula  [26]  and  recalling  Equation 


i 
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(3.68),  (3.66)  and  (3.65),  the  six  roots  for  Equation  (3.63)  are 
found  to  be 


Ao(l,2)  =  1  1.2364  ao 
Xo(3,4)  =  ±  0.2903  ao 
Xo(5,6)  =  ±  0.9659  Uq 


(3.71) 


Hence  the  zeroth  order  composite  solution  (Section  1.3)  takes 
the  form 

fkl 

Uc^^'  =  ((Pie  ’  °  +  Pge  2  °  +  Pae  ^  )cos  aoZ)  ^ 

=  -  0.7433  ExbY  +  {(Rie  '  °  +826^^° 


,(•<) 


+  830  )C0S  aoZ} 

(k)  -6o««n 

\^/  _  rk  rsrrkA  i  r  f  n  ^  ^ 


(3.72) 


0.0604  exhZ  +  {($16  "  +$26 

-8,a  n  (1(1 

+  $30  ^  °  )sin  aoZ} 


7“o 


(k) 


where  Bi  =  1.2364 
62^*^^  =  0.2903 
83^*^^  =  0.9659 

Satisfying  the  governing  equations  and  the  boundary  conditions 


leads  to  the  following  equations: 
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Now  consider  Equations  (3.69).  It  is  clear  that  Layer  1 
(+  45°)  and  Layer  2  (-  45°)  are  antisyimetric  in  U  and  symmetric  in 
V  and  W  with  respect  to  the  infinitesimal  thin  slice  (Fig.  17d). 
Upon  enforcing  exact  continuity  in  displacements  at  Z  =  »  the 


lim  cos 

C-K) 


1 

a 

0  .  (1). 

=  lim  cos 

°  -  a  <2). 

2  ^ 

2  “o'* 

C-K) 

=  0 


(3.76) 


which  gives 


=  (2n  +  l)'!i  t  n  -  0,1»2,.... 
=  (2n  +  l/iT  ,  n  =  0,l»2>.... 


(3.77) 


Hence , 


cos 


(1) 

(2)  1 

“0  +  r  (1) 

=  COS 

"0  -  c  a 

2  “o 

2  ^  “o 

J 

-  0 


(3.78) 


for  0  <  c  «<  1 


where  and  are  given  in  Equations  (3.77), 

f.  W 


Thereby 


lim 

C-H) 


tan 


(3.79) 


and 


tan 


L  (k) 


2  -“o'^ 


=  K 


(3.80) 


0  <  ^  <<<  1 


where  K  is  a  finite  large  posicive  value. 
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At  this  stage,  assigning  any  large  value  for  K  determines  the 


corresponding 


and  hence  the  interlaminar 


stresses.  It  may  be  shown  that  the  only  stresses  related  to 


tan 


'a„(k) 


±  a, 


(k), 


(k) 


are  and 


(k) 


However,  the  latter 


2  ^  ^0 

anishes  identically  at  the  free  edge  as  required  by  the  stress  free 
oundary  conditions  (2.17).  Hence  the  singular  behavior  is  found  in 
at  the  intersection  of  the  free  edge  and  the  interfacial  plane 
=  -  .  This  provides  a  definite  mathematical  evidence  for  the  pre- 
icted  singularity  in  Reference  [7]  and  will  be  further  discussed  in 

he  following  chapter. 

The  interlaminar  stressea  re  plotted  in  Figures  18  -  20. 

.2.2  [-45/4515  GRAPHITE-EPOXY  LAMINATE 

Consider  the  laminate  of  Fig.  17e.  Interchanging  the  super- 
cripts  1  and  2  in  Equations  (3.62)  through  (3.73)  gives  a  composite 
lolution  identical  to  (3.72).  The  corresponding  interlaminar  stresses 


are  shown  in  Figures  21  ~  23. 


NTERLAMINAR  shear  stress  Tyz  IN  [45/-452 
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FIGURE  21.  INTERLAMINAR  SHEAR  STRESS 
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INTERLAMINAR  NORMAL  STRESS  CT^  IN  [-45/45] 


Chapter  IV 

RESULTS  AND  DISCUSSION 

In  t!.e  preceding  chapter,  the  general  method  of  solution  of 
Chapter  II  was  applied  to  the  special  graphite-epoxy  laminates  [O/DOlj, 
[SO/Olj,  [45/-45]g  and  [-45/45jg.  To  demonstrate  the  capability  of 
the  solution  the  results  for  these  laminates  are  presented  and  dis¬ 
cussed  in  this  chapter. 

4.1  THE  FOUR  LAYER  UNIDIRECTIONAL  LAMINATES 

It  has  been  stated  in  Section  1.3  that  the  accuracy  of  the  pertur¬ 
bation  solution  depends  upon  the  perturbation  parameter  c.  That  is, 
the  smaller  e,  the  better  the  result.  This  will  be  demonstrated  in 
what  follows. 

The  interlaminar  shear  stress  and  the  interlaminar  normal 
stress  a2  (the  peel  stress)  as  functions  of  the  perturbation  parameter 
e  are  presented,  respectively,  in  Figures  8  and  9  for  the  [f/30]^ 
laminate.  From  the  figures,  it  is  clear  that  the  boundary  layer  width 
becomes  smaller  as  e  decreases  in  magnitude.  (Asymptotic  recovery  of 
the  lamination  theory  is  implied  by  the  incomplete  domain  of  ^.1  It 

D  ■ 

should  be  noted  that  the  relative  extreme  values  of  the  stresses  are 
finite  and  remain  unchanged  as  c  decreases.  This  indicates  that  the 
present  theory  is  capable  of  approximating  the  maximum  value  of  the 
interlaminar  stress  intensities  for  intermediate  as  well  as  small 
values  of  e.  Also,  the  difference  betv/een  the  cases  c  =  0.133  and 
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G  =  0.050  is  much  more  than  that  between  g  =  0.050  and  g  =  0.033. 

Vihile  the  curve  of  g  =  0.033  serves  as  the  most  accurate  of  the  three 
stress  results  for  their  corresponding  g,  it  lends  confidence  to  say 
that  for  this  [0/90]s  graphite-epoxy  laminate,  a  geometric  ratio  of 
0.050  (=  Iq-)  or  smaller  is  sufficiently  small  to  lead  to  good  results 
using  the  present  method  of  solution. 

Numerical  results  obtained  by  this  author  using  the  finite  dif¬ 
ference  program  of  Pipes  [12]  indicate  that  the  smallest  geometric 
ratio  for  which  that  program  gives  physically  admissible  result  is 

c  =  0.133  (*77).  Below  this  ratio,  the  instability  in  the  solution 
1 5 

does  not  yield  satisfaction  of  the  force  equilibrium  Fy  =  0  (Fig.  2 
and  Equation  (2.26)).  This  may  be  attributed  to  the  inherent  sensi¬ 
tivity  of  the  finite  difference  approximation  to  the  ratio  of  the  grid 
spacings  for  partial  differential  equations  [27], 

Comparisons  between  the  results  of  the  finite  difference  solution 
and  the  present  theory  are  presented  in  Figures  10  and  11  for  the  case 
G  =  0.133.  From  Figure  10  it  is  clear  that  the  present  theory  tends  to 
predict  a  higher  maximum  intensity  for  the  interlaminar  shear  stress 
Ty^.  The  boundary  layer  width  is  approximately  the  same  for  both 
solutions.  Figure  11  shows  that  the  present  theory  predicts  a  smooth, 
continuous  distribution  for  02  which  identically  satisfies  the  self- 
equilibrating  condition  l  F^  =  0  (Equation  (2.70))  whereas  the  finite 
difference  solution  yields  unstable  results  near  the  free  edge  which 
obviously  do  not  satisfy  this  equilibrium  condition.  In  regions 
removed  from  the  free  edge,  both  solutions  indicate  asymptotic  recovery 


of  the  lamination  theory. 

Figures  12  and  13  show  the  interlaminar  stresses  for  the  [QO/Ojj 
laminate"the  reversed  stacking  sequence  from  the  previous  example. 

From  the  figures,  the  physical  validity  of  the  present  theory  is  con¬ 
firmed  by  the  sign  reversals  in  both  Ty^  and  as  a  result  of  force 
and  moment  equilibrium  (Fig.  7).  Again,  boundary  layer  shifts  due  to 
the  reduction  in  e  are  observed.  The  maximum  stress  intensities  of 
Ty2  and  oj.  in  the  C90/0]5  laminate  are  found  to  be  finite  but  higher 
than  those  in  the  [O/QO]^  laminate  (Figs.  8  and  9).  This  is  due  to 
the  fact  that  in  the  calculation  for  the  [QO/Ojg  laminate,  the  0°  layer 
was  employed  as  the  reference  layer.  On  the  other  hand,  in  the  calcu¬ 
lation  for  the  [0/90]^  laminate,  the  90°  layer  was  employed  as  the 
reference  layer. 

Comparisons  between  the  finite  difference  results  and  the  present 
theory  are  presented  in  Figures  14  and  15.  The  present  theory  again 
predicts  a  higher  Ty^  than  the  finite  difference  solution.  Also,  the 
present  theory  yields  a  more  acceptable  distribution  for  the  inter¬ 
laminar  normal  stress  in  view  of  the  zero  stress  resultant  require¬ 
ment.  In  regions  removed  from  the  free  edge,  the  lamination  theory  is 
recovered  asymptotically  in  both  solutions. 

4.2  THE  FOUR  LAYER  ANGLE-PLY  LAMINATES 

Pipes  and  Pagano  [7]  pointed  out  that  the  interlaminar  shear 
stress  in  a  [ASZ-ASjg  laminate  tends  to  grow  without  bound  near 
the, free  edge  (Section  1.1).  Hence  the  calculated  maximum  intensity  of 
Ty2  by  the  finite  difference  approximation,  though  higher  than  other 
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numerical  investigations  [6,  8.  11],  is  still  very  questionable.  It 
was  discussed  in  Section  1.2  that  failure  to  satisfy  some  stress  free 
boundary  conditions  were  observed  in  the  finite  difference  solution. 
Also,  these  results  showed  no  sign  reversals  in  oy.  and  ty^  in 
consequence  of  reversing  the  stacking  sequence. 

In  the  present  theory  the  mathematical  evidence  for  the  singu¬ 
larity  in  Txz  can  be  shown  (Subsection  3.2.1)  to  be  in  terms  of 
the  trignometric  equation 

(2) 

tan  (^^2  -  =  K 

where  0  <  c  «<  l  and  K  is  a  near-singular  large  number.  The  value  of 
(2) 

^  must  satisfy  equation  (3.75) 

tan  (2^  . 

jYy - "  0-5  ,  0  <  c  «<  1 

(4.2) 

Obviously,  the  limiting  analysis  of  the  present  theory  (Fig.  ]7d,  h) 
provides  no  unique  determination  of  the  value  of  K.  It  is  only  through 
experimental  investigation  that  this  value  may  be  realistically  deter¬ 
mined.  Such  an  investigation  should  be  considered  as  a  future  study. 
For  the  purpose  of  comparisons,  K  is  taken  to  be  20.3713,  a  value  that 
leads  to  a  maximum  stress  intensity  within  the  elastic  limit. 

Comparisons  between  the  results  of  the  finite  difference  solution 
and  the  present  theory  are  presented  in  Figures  18  through  23.  Figure 
18  shows  the  variation  of  the  interlaminar  shear  stress  the 


interface  2  =  ^  .  The  near-singular  free  edge  intensity  of  the  present 
theory  is  much  higher  than  the  finite  difference  result  and  the 
boundary  layer  width  is  much  smaller.  Figure  19  shows  the  variation 
0^  the  interlaminar  shear  stress  Both  solutions  satisfy  the 

stress  free  boundary  condition  Xy^  =  0  at  the  free  edge.  The  negative¬ 
positive  variation  of  the  present  theory  confirms  the  additional  self- 
equilibrating  condition 

^  fy  ~  f  Xy2  dy  =0  (4.3) 

■'o 

(as  a  result  of  the  zeroth  order  vanishing  of  oy  in  the  interior 
region).  The  finite  difference  solution,  on  the  other  hand,  cannot 
satisfy  such  a  condition.  The  erroneous  of  the  finite  dif¬ 
ference  solution  at  the  free  edge  (not  shown  in  the  figures),  as 
described  in  Section  1.2,  is  believed  to  be  caused  by  inherent  errors. 

In  Figure  20  the  interlaminar  normal  stress  of  the  finite  difference 
solution  indicates  instability  near  the  free  edge;  hence,  no  comparison 
can  be  made  betvieen  the  two  solutions  in  this  region.  Since  the  auto¬ 
matic  satisfaction  of  the  self-equilibrating  condition 

I  F^  =  dy  =  0  (4.4) 

•'o 

has  been  demonstrated  by  the  present  theory  (Chapter  III)  and  can  be 
observed  from  the  figure,  and  since  is  not  proportional*  to  the 


*As  shown  in  Equations  (3.46)  -  (3.53). 
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near-singular  value  of  K,  the  present  theory  is  believed  to  have  pre¬ 
dicted  a  more  accurate  maximum  finite  intensity  of  the  interlaminar 
normal  stress.  Such  a  determination  is  most  important  in  the  delami¬ 
nation  failure  mode  [21,  28]  of  composites.  Although  the  moment  self- 
equilibrating  condition  (Fig.  16)  is  not  directly  observable  from 
Figure  20,  the  magnitude  of  this  couple  moment  can  be  determined  as 

“  =  y  dy  -  0.0027  ^5y2[liK5ff)  W-5) 

f  2 ) 

where  a  equals  2K,  a  near-singular  value  from  Equations  (4.1)  and 
(4.2).  Hence  the  self-equilibrating  CL..dition  of  the  couple  moment  is 
confirmed  immediately. 

When  the  stacking  sequence  is  reversed  to  [-AS/AS]^  (Figs.  21,  22, 
23),  the  interlaminar  shear  stress  experiences  a  sign  change  in 
order  to  balance  the  central  plane  shear  resultant  Txy^^^(0,Z)hdZ 

2 

which  also  experiences  a  sign  reversal.  The  sign  of  both 
remain  unchanged.  This  is  in  agreement  with  the  finite  difference 
results  (Figs.  19  and  22,  20and  23).  For  e  =  0.133,  the  finite  dif¬ 
ference  solution  predicts  a  small  uniform  oy  along  the  central  plane 
(not  shown  in  the  figures)  which  does  not  change  its  sign  and  magnitude 
for  the  reversed  stacking  sequence.  For  c  £0.0133  the  finite  dif¬ 
ference  solution  yields  erroneous  results  for  oy  due  to  the  instability 
of  the  solution.  The  present  theory  exhibits  no  such  instabilities. 

It  is  important  to  note  that  the  interlaminar  normal  stress 
is  independent  of  the  stacking  sequence  and  always  tensile  near  the 
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free  edge.  For  both  the  [45/-45]5  and  the  [-45/45]5  laminates,  a 
finite  maximum  intensity  is  predicted  at  the  exact  free  edge  (Figs.  20, 
23).  This  indicates  that  the  oelamination  failure  mode  [21,  28]  should 
always  be  considered  for  reliable  design  of  such  laminate  configura¬ 
tions. 

It  is  clear  that  the  present  study  has  obtained  improved  results 
for  the  interlaminar  behavior  of  the  [45/-45]5  and  [-45/45]5  graphite- 
epoxy  laminates.  Since  the  aforementioned  self-equilibrating  conditions 
were  originally  considered  for  the  2m  layer  angle-ply  laminate  (Fig. 

16),  the  interlaminar  stress  variations  in  any  angle-ply  laminate  may 
be  expected  to  be  similar  to  those  in  Figures  18  through  23. 

4.3  ACCURACY  AND  LIMITATIONS 

As  discussed  earlier,  the  accuracy  of  the  present  theory  depends 
upon  the  geometric  ratio  c  =  ^  .  Hence,  the  relative  order  of  magni¬ 
tude  of  the  individual  terms  in  the  governing  equations,  in  relation  to 
,  should  be  further  discussed. 

4.3.1  BIDIRECTIONAL  LAMINATES 

The  coupled  governing  differential  equations  for  bidirectional 
laminates  (Equations  (3.2))  are 

(Q22t5')^V,yy  +  Q44V,22  +  (Q44  +  Q23)(^)W.yz  = 

(4.6) 

{(Q44  +  Q23)(F^V,y2  +  Q44(^)^W,yy  +  0331^,2^  =  0}^''^ 

From  these  dimensionless  equations,  it  is  essential  that  the 
order  of  magnitude  of  the  coefficients  of  V,yy  and  W,yy,  V,y2  and 
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M,yz»  V.zz  and  be  0(e  )  <  0(F)  <  0(1)*,  respectively,  in  order  to 
properly  stretch  the  boundary  layer  region  with  a  transfomiation  in  the 


(4.7) 


Hence,  if  the  material  properties  are  fixed,  the  geometric  ratio  r 

b 

obviously  plays  the  dominant  role.  For  the  graphite-epoxy  laminate 
with  g- =  0.133  (Chapter  III),  Equations  (4.6)  may  be  transformed  to 


(90®) 


0.046  V,yY  +  V,zz  +  0.208  W,y2  =  0 
0.080  V,Yz  0.007  Wiyy  ^  ^’ZZ  ~  ^ 


Jo. 42  V,YY  +  V,zz  +  0.22  W,^ 


0.085  V,Yz  0.007  Wjyy  W,zz  ”  0 


It  may  be  observed  that,  for  this  geometric  ratio,  a  perturba¬ 
tion  solution  using  the  0®  layer  as  the  reference  layer  should  lead  to 
more  accurate  results. 

If  the  geometric  ratio  is  now  reduced  to  0.050  for  the  same 
laminate  material,  Equations  (4.6)  become 


(90®) 


0.0065  V,yy  +  V.zz  +  0.0782  W.y^  = 
0.030  V,y2  +  0.001  W,yy  +  W,zz  ”  0 

0.05  V,yy  +  V,zz  0.083  Wjyz  ~  0 

0.032  V,yz  +  0.001  W.yy  +  W,2z  '  0 


(4.10) 


(4.11) 


*  G  represents  the  approximate  order  of  the  products  of  0- 

and  G. 
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It  is  clear  that  the  order  of  each  term  relative  to  ^,22  or 
W,22  shrinks  as  e  diminishes.  This  means  that  the  degree  of  accuracy 
of  the  zeroth  order  perturbation  solution  is  improved  by  the  diminish¬ 
ing  geometric  ratio.  Obviously,  these  equations  provide  mathematical 
evidence  to  support  the  judicious  statement  made  in  Section  4.1— 

"A  geometric  ratio  of  0.050  or  smaller  leads  to  sufficiently  accurate 
results." 

4.3.2  ANGLE-PLY  LAMINATES 

The  governing  differential  equations  for  the  i  45°  graphite-epoxy 
laminate  with  0.133  are 

0-111  U.yy  +  U,22  -  0.094  V,yy  -  0.007  W.y^  =  0 
-0.094  U,yy  +  0.140  V,yy  +  V,22  +  0,215  W,y2  =  0  (4.12) 

-0.003  U,y2  +  0.083  V,y2  +  0.007  W,yy  +  W,22  =  0 

O.m  U.yy  +  0,22  0.094  V,yy  +  0.007  W,Y2  =  0 

(-45°;  0.094  U,yy  0.140  V,yy  +  V,22  +  0.215  W.yz  =  0  (4.13) 

0.003  U,Yz  +  0.083  V,y2  +  0.007  W.yy  +  W,2Z  =  0 

Again,  the  order  of  magnitude  of  each  coefficient  relative  to 

W,22  can  be  observed.  From  Equations  (4.12)  and  (4.13), 
it  can  be  safely  stated  that  the  present  theory  should  lead  to  suf¬ 
ficiently  accurate  results  for  =  0.133  or  below.  Hence  no  further  ' 
reduction  of  the  geometric  ratio  needs  be  elaborated  upon. 

From  the  above  discussion,  the  fact  that  the  boundary  layer 
penetration  becomes  weaker  as  ^  decreases  can  be  detected  simply  by 
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examining  the  relative  orders  of  various  terms  in  the  governing  dif¬ 
ferential  equations.  The  smaller  the  coefficients  relative  to  the  co¬ 
efficient  one  of  U»2Z»  ^^21  *^’ZZ»  weaker  the  boundary  layer 

effect.  Hence  the  uniform  stress  distribution  in  the  central  plane 
(Section  2.2)  is  justified  by  the  more  rapid  recovery  of  the  lamina¬ 
tion  solution. 

It  must  be  recalled  that  in  the  interior  region  of  the  present 
theory  the  exact  satisfaction  of  the  vanishing  stress  boundary  condi¬ 
tions  on  the  top  and  bottom  surfaces,  the  continuities  in  the  inter¬ 
laminar  stresses,  and  the  force  equilibrium  in  the  central  plane  were 
inforced.  Also  in  the  interior  region  the  exact  displacement  continuity 
in  U  and  V  were  satisfied  by  the  modified  zeroth  order  interior 
region  (Subsection  2.3.2).  For  bidirectional  laminates,  the  slight  dif¬ 
ference  found  in  the  displacement  W  (Table  1)  for  the  two  layers  may 
be  reduced  or  eliminated  by  higher  order  considerations.  This  is 
mainly  due  to  the  differences  in  material  properties  that  constitute 
the  governing  differential  equations.  For  the  angle-ply  laminates, 
the  exact  continuity  in  this  displacement  v;as  found  to  be  satisfied 
automatically  (Table  2). 

In  the  boundary  layer  region,  the  bidirectional  laminates  again 
reveal  differences  in  the  exponential  functions  (Equations  (3-33)) 
owing  to  the  intractable  material  dissimilarities.  Hence  no  exact 
displacement  continuity  in  this  region  may  be  inforced  for  the  limit¬ 
ing  free  body  considered  in  Figure  7.  Nevertheless,  the  satisfaction 
of  the  symmetry  conditions,  the  stress  boundary  conditions  at  fie  free 


edge,  the  force  and  moment  equilibrium  about  this  free  body,  and  the 
continuity  in  interlaminar  tractions,  is  believed  to  have  attained  a 
solution  which  is  an  improvement  over  previously  available  solutions. 

On  the  other  hand,  the  boundary  layer  region  solution  for  the 
angle-ply  laminate  exactly  satisfies  the  symietry  conditions,  the 
stress  boundary  conditions  at  the  free  edge,  the  displacement  and 
stress  continuity  conditions,  and  the  force  and  moment  equilibrium  re¬ 
quirements  (Figs.  16,  17).  This  exactness  of  the  present  theory  is 
entirely  due  to  the  favorable  parametric  relations 


C  (1)  -  r  (2) 

^ij  ~  ^ij  » 

i  =  1,  2,  3  and  j  =  1,  2,  3 

Ckk'*  = 

k  =  4,  5,  6 

(4.14) 

r  (1)  .  .  (2) 

^n6  ~  “^n6  » 

n  =  1,  2,  3 

Moreover,  it  is  this  exactness  that  leads  to  the  mathematical  evidence 
for  the  stress  singularity  in  the  interlaminar  shear  stress 

4.4  GENERALITY  AND  APPLICABILITY 

From  the  solution  method  developed  in  Chapters  II  and  III,  it  is 
clear  that  the  detailed  solution  procedures  of  the  present  theory  can 
be  readily  programmed  for  a  computer.  The  simple  r.alculatiori  steps 
require  no  approximate  or  iterative  techniques.  The  generality  of  the. 
theory  can  be  directly  applied  to  variable  laminate  configurations 
(Fig.  1)  with  more  layers  than  any  existing  computer  program  can 
possibly  handle,  thermal  strains  can  be  readily  included  through  the 
constitutive  equations  (2.1)  to  determine  the  induced  thermal  stresses 
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due  to  the  temperature  drop  from  the  curing  temperature  of  the  laminate 
or  for  laminater,  in  a  high  temperature  environment.  The  computer 
program  thus  developed  would  be  capable  of  predicting  the  interlaminar 
stress  intensities  between  any  two  layers  including  the  midplane  Z  =  0. 

» 

Failure  hypotheses  can  then  be  established  based  on  the  interfacial 
plane  where  the  stress  intensity  reaches  a  relative  maximum. 

Finally  the  generality  of  the  solution  procedures  in  the  present 
theory  can  be  directly  applied  to  explore  related  problems  such  as  a 
laminate  with  internal  free  edge  in  the  form  of  center  holes,  cracks, 
etc.,  a  laminate  subject  to  pure  bending  at  the  ends  x  =  ±  L,  time- 
depondent  boundary  layer  effects  due  to  cyclic  loadings,  and  so  forth. 
The  important  experimental  determination  of  the  material  parameter  K 
may  also  be  pursued  as  a  future  research  topic. 
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Chapter  V 
CONCLUSIONS 

In  the  present  thesis  a  general  method  of  solution  for  a  balanced 
symmetric  composite  laminate  subject  to  a  uniaxiai  extension  has  been 
developed  based  upon  a  perturbation  analysis  of  an  elastic  limiting 
free  body  containing  an  interfacial  plane. 

In  summary  of  the  theoretical  achievements  of  the  present  study, 

the  following  conclusions  can  be  made. 

(1)  The  solution  satisfies  the  symmetry  conditions,  the  stress 
free  boundary  conditions,  most  continuity  conditions,  and 
the  force  and  moment  equilibrium  of  the  limiting  free  b:dy. 

(2)  The  solution  predicts  smooth  continuous  interlaminar 
stresses  with  no  instabilities. 

(3)  The  solution  provides  the  finite  maximum  magnitude  of  the 
interlaminar  normal  stress  for  all  laminate  configura¬ 
tions. 

(4)  For  given  material  properties,  the  solution  accuracy  depends 
upon  the  geometric  ratio  e  .  For  [0/90]s  Gr/E  laminate. 

'  <  I_  leads  to  satisfactory  results  while  for  C±453s  Pr/E 

-  20 

laminate,  e  predicts  satisfactory  results. 

(5)  For  all  laminates  with  geometric  ratio,  0  <  g  «  1,  high 
gradient  displacement,  strain  and  stress  fields  are  shown 
to  exist  near  the  free  edge. 


(6)  The  above  boundary  layer  effect  decays  exponentially  to 
recover  the  lamination  solution  in  the  interior  regions. 

(7)  For  bidirectional  laminates  the  axial  displacement  function 
U  is  identically  zero.  Hence  no  or  xxz  exist  in  the 
laminate. 

(8)  For  bidirectional  laminates,  the  interlaminar  normal  stress 
©2  is  finite  with  the  sign  depending  upon  the  stacking 
sequence.  For  example,  for  a  [O/QO]^  Gr/E  laminate,  a 
maximum  tensile  exists  at  the  free  edge  v/hile  for  a 
[GO/Ojj  Gr/E  laminate,  a  maximum  compressive  is  predicted. 

(9)  For  angle-ply  laminates,  the  exactness  of  the  solution  leads 
to  the  mathematical  evidence  of  singular  interlaminar  shear 
stresses  and  Xy^  at  or  near  the  free  edge. 

(10)  For  angle-ply  laminates,  the  interlaminar  normal  stress 
takes  on  a  finite  maximum  tensile  value  at  the  free  edge,  and 
is  independent  of  the  stacking  sequence. 

(11)  The  solution  procedure  can  be  readily  prograrrmed  for  a 
computer.  Such  a  generalized  computer  program  would  be 
capable  of  predicting  interlaminar  stresses  between  any 
two  layers  of  a  general  multi-layered  laminate. 

(12)  The  present  theory  suggests  vaiable  means  for  solving 
important  related  problems  of  practical  interest. 
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